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PREFACE 


The object of this book is to provide a compact ex¬ 
position of the fundamental resuits in tiie theory of 
tensors and also to illustrate the power of the tensor 
technique by applications to differential geometry, 
elasticity, and relativity. In the first five chapters the 
mathematical concepts are developed without undue 
stress on rigour. The remaining three chapters are in¬ 
dependent of one another except that sections 88 and 89 
of chapter VI, which treats Euclidean three-dimensional 
differential geometry, arc necessary for a proper under¬ 
standing of chapter VII which contains the theory of 
cartesian tensors and elasticity. Finally, chapter VIII is 
devoted both to the special and general theories of 
relativity* In the limited space available it is impossible 
to do justice to the physical principles underlying both 
these theories* But in order to help the reader unac¬ 
quainted with relativity some explanatory matter has 
been incorporated into the text. 

The presentation owes much to the authors listed in the 
bibliography, especially to McConnell, Synge and Schild. 
In particular, I wish to express my thanks to Dr, D. E* 
Rutherford for numerous suggestions and helpful 
criticisms during the manuscript and proof stages* Lastly 
I wish to thank my wife for her help with the proof¬ 
reading. 

D* S* 

Trinity College, Dublin; July, 1952 

PREFACE TO THE THIRD EDITION 

In this edition various errors have been corrected. 
Further, I wish to thank Mr. L. Lovitch for a neat proof 
{inserted on page 56) that the curvature tensor is zero in 
a flat space* 


B* S. 
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CHAPTER I 


TENSOR ALGEBRA 

§ 1. Introduction 

The concept of a tensor has its origin in the develop* 
ments of differential geometry by Gauss, Riemann and 
Christoffel. The emergence of Tensor Calculus, otherwise 
known as the Absolute Differential Calculus, as a syste¬ 
matic branch of Mathematics is due to Ricci and his pupil 
Levi-Civita. In collaboration they published the first 
memoir on this subject: — 'Methodes de calcul differential 
absolu et leurs applications* y Maihematische Annalen, vol. 
54, (1901). 

The investigation of relations which remain valid when 
we change from one coordinate system to any other, is 
the chief aim of Tensor Calculus. The laws of Physics 
can not depend on the frame of reference which the phy¬ 
sicist chooses for the purpose of description. Accordingly it 
is aesthetically desirable and often convenient to utilise 
the Tensor Calculus as the mathematical background in 
which such laws can be formulated. In particular, 
Einstein found it an excellent tool for the presentation of 
his General Relativity theory. As a result the Tensor 
Calculus came into great prominence and is now invalu¬ 
able in its applications to most branches of Theoretical 
Physics; it is also indispensable in the differential geome¬ 
try of hyperspace. 

It is assumed that the reader has an elementary know¬ 
ledge of determinants and matrices. As he may not be 
acquainted with the Calculus of Variations, the minimum 
problem in the theory of geodesies is treated from first 
principles. 
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§ 2, N-Dimensional space 

Consider an ordered set of N real variables x 2 t . , 

. .,, xF; these variables will be called the coordinates 
of a point, (The suffixes I, 2, * , * i, * * ■ N, which we shall 
call superscripts, merely serve as labels and do not possess 
any significance as power indices. Later we shall introduce 
quantities of the type a £ and again the i , which we shall 
call a subscript, will act only as a label,) Then all the 
points corresponding to all values of the coordinates are 
said to form an N-dimensional space, denoted by Vn* 
Several or all of the coordinates may be restricted in 
range to ensure a one-one correspondence between points 
of the V N and sets of coordinates, 

A curve in the V N is defined as the assemblage of 
points which satisfy the N equations 

x € — (f = 1, % . *, N) 

where u is a parameter and x i ( , u) are N functions of u, 
which obey certain continuity conditions. In general, it 
will be sufficient that derivatives exist up to any order 
required, 

A subspace V M of V N is defined for M < A 7 as the 
collection of points which satisfy the AT equations 
x i ss x^u 1 , u 2 , * , . u M ), (i = I, 2, ,, , A T ) 
where there are M parameters u l , u 2 , *,, u M , The 
* , , u M ) are N functions of the u 1 , u 2 , * , * u M 
satisfying certain conditions of continuity. In addition 
the M X N matrix formed from the partial derivatives 
dx i jdvp is assumed to be of rank M *, When M = A T — 1* 
the subspace is called a hypcrsurface. 

§ 3, Transformation of coordinates 

Let us consider a space Uv with the coordinate system 
X 1 , x 2 , , , , The N equations 
(3,1) «* = x 2 , . flfW), (i = 1, 2, ,., AT) 

* T, Levi-Ctvita, The Absolute Differential Calculus, pp, 9—12, 
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§4 

where the q> £ are single-valued continuous differentiable 
functions of the coordinates, define a new coordinate 
system x l 7 x 2 , , . , 0. Equations (3,1) are said to define 
a transformation of coordinates. It is essential that 
the N functions be independent, A necessary and suf¬ 
ficient condition is that the Jacobian determinant formed 
from the partial derivatives dx i jdx i does not vanish*. 
Under this condition wc can solve equations (3*1) for the 
x* as functions of the x 1 and obtain 

x* — ^(x 1 , x\ ... x N ) (i = 1, 2,... N). 

§ 4, Indicial and summation conventions 
We will now introduce the following two conventions: 

(1) , Latin indices , used either os subscripts or superscripts, 
will take all values from l to N unless the contrary is 
specified. Thus equations (3,1) are briefly written 

— <p i (x\ x\ , ,, x^) ? the convention informing us that 
there are N equations, 

(2) , Jf a Latin index is repeated in a term , then it is un¬ 
derstood that a summation with respect to that index over 
the range 1 , 2, ,, , N is implied , Thus instead of the cx- 

N 

pression Sa i x i % we merely write a i x\ 

i—1 

Now differentiation of (3,1) yields 
3 CD* at 3~i 

which simplify, when the above conventions are used, to 

(^*1) dx* — dx r , 

dx T 

I he repeated index r is called a dummy index, as it 
can be replaced by any other Latin index, except i in 

* R. P, Gillespie, Partial Differentiation, pp, 43—4G, 
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§ 5 

this particular case. That is, equations (4.1) can equally 

30 

well be written d& = ^ dx m or for that matter 

pi-j 0X m 

A 0 = _ d#*. In order to avoid confusion, the same 
dx i % , 

index must not be used more than twice in any single 
N , , . 

term. For example (2J a^x *) 2 will not be written a^x d { x 

i^i 

hut rather a i a^x i xK It will always be clear from the con* 
text whether x 2 means x with superscript 2 or x squared. 
Usually powers will be Indicated by the use of brackets; 
thus (a ^) 2 means the square of scP. The reason for using 
superscripts and subscripts will be indicated in due 
course. 

Let us introduce the Kronecker delta defined by 

r df = i if # = * 

( 4 - 2 ) l«J- o if j=£Je. 

An obvious property of the Kronecker delta is that 
SfA 1 = A*, since in tiic left-hand side of this equation 
the only surviving term is that for which j = k. Also 
3 x k jdx* = 6 f, because the coordinates x i are independent, 

Ex. Show that 6 j$ k = <S*; = N; — = °i' 


§ 5. Contravariant vectors 

A set of N functions A 1 of the N coordinates x l are said 
to be the components of a contravariant vector if they 
transform according to the equation. 


(5.1) 


A< = f-, AK 
dx j 


on change of the coordinates x* to x** Ibis means that 
any N functions can be chosen as the components of a 
contravariant vector in the coordinate system x i , and the 
equations (5*1) define the N components in the new 
coordinate system x*. On multiplying equations (5.1) by 
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dx^jdx 1 and summing over the index i from 1 to A 7 , we obtain 
dx k dx* dx i dx k k 


Hence the solution of equations (5*1) is 
(5.2) 


.. dx k T 

A ' " 5? *■ 


When we examine equations (4.1) we see that the 
differentials dx i form the components of a contravariant 
vector, whose components in any other system are the 
differentials dx* of that system* It follows immediately 
that dx i jdu is also a contravariant vector, called the 
tangent vector to the curve x* — ^(u). 

Consider now a further change of coordinates 
x H = g l (x\ z 2 ,... x N ), Then the new components A n 
must be given by 

AH _ 3®' 1 jj _ dx’i 30 _ dx H 

A ~ dx* A ~ 30 dz k A dx* ‘ 


This equation is of the same form as (5.1), which shows 
that the transformations of contravariant vectors form 
a group ** 

With the exception of the coordinates x i themselves 
a single superscript will always denote a contravariant 
vector unless the contrary is explicitly stated. The coor¬ 
dinates x l will only behave like the components of a 
contravariant vector with respect to linear transfor¬ 
mations of the type x i — a]x^ where the aj are a set 
of A' 2 constants, which do not necessarily form the com¬ 
ponents of the entity introduced in section 8 and there 
called a tensor. For in this ease dx i jdx i — a) and the trans- 

formation can be rewritten — -x—- xK W r itli respect to 

ox* 

general transformations of coordinates, the x l do not form 


* W. Ledennann, Introduction to the Theory of Finite Groups, 
PP- 2^3. 
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the components of a contravariant vector. Essentially 
this means that if we select A i = x\ then the new com¬ 
ponents A i with respect to the coordinate system 0 do 
not satisfy the equations A* — x i * 


Ex, If a vector has components 


d*x d*y 
dF’ dt* 


in rectangular cartesian 


dH IdOy d 2 Q 2 dr dO . 

coordinates, show that they arc — -r ^ j , — + - J t J ( »» 
polar coordinates. 


§ 6. Covariant vectors 

A set of N functions of the N coordinates x' are said 
to be the components of a covariant vector if they 
transform according to the equation 

(6-1) 


on change of the coordinates 0 to 0. Any N fimetions 
can be chosen as the components of a covariant vector 
in the coordinate system x i , and the equations {6,1) 
define the N components in the new coordinate system x\ 
On multiplying (6,1) by d&jdx* and summing over the 
Index i from 1 to A T , we obtain 


( 6 . 2 ) 


d0 j d0 dx* _ dx* . _ , 
dx k 1 dx k d0 i dx k * k 


Since — — it follows immediately from (6,1) 

ox 1 ox ! 

that the quantities dfjdx 1 are the components of a co¬ 
variant vector, whose components in any other system 
are the corresponding partial derivatives 3// d0. Such a 
covariant vector is called the gradient of /. 

A single subscript will always denote a covariant vector 
unless the contrary is explicitly stated. In conformity 
with this convention we shall regard the index i in the 
covariant vector djjdx 1 as a subscript. 
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We now show that there is no distinction between con* 
travariant and co variant vectors when we restrict our¬ 
selves to transformations of the type 

(6.3) 0 ~ a f m x m + b 

where b £ arc A T constants which do not necessarily form 
the components of a contravariant vector and a l m are 
constants (not necessarily forming a tensor) such that 

oUL = > 

r m m 

We multiply equations (6,3) by a* and sum over the 
index i from 1 to N and obtain 


Thus 


x r = a*0 — 
dx l dx i 


which shows that the equations (5,1) and (6,1) define 
the same type of entity. 

Ex, Prove that the transformations of co variant vectors form a 
group. 


§ 7. Invariants 

Any function I of the N coordinates x l is called an 
invariant or a scalar with respect to coordinate trans¬ 
formations if I — /, where / is the value of / in the new 
coordinate system 0 . 

I rom the components A i and of a contravariant and 
covariant vector respectively, we can form the sum A ( B ir 
When we change to new coordinates 0 t this sum trans¬ 
forms to AtB^ Now 


IiB( = S^S^= q Ai b *= AkB *■ 

That is, 

Am, = A*B t . 

Thus A t B i is an invariant. 
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Another invariant is the quantity 

& * = ai + aS + ... + a£«ff. 


§ 8. Second order tensors 

Form the A 72 quantities A fi — B'CJ, where B* and C* 
are tile components of two contravariant vectors. It 
follows from (5.1) that the transform according to 


( 8 . 1 ) 


, dx 1 dx> 

'**-©» A - 


More generally, if we have N 2 functions A if whose 
transformation law is that of (8.1), then we call A w the 
components of a contravariant tensor of the second 
order. This tensor is not necessarily the product of the 
components of two contravariant vectors. Any set of A 1 
functions can be chosen as the components of a contra¬ 
variant tensor of the second order, and then (8,1) defines 
the components in any other coordinate system 0. 

Similarly* if we have N 2 functions A iS whose transfor¬ 
mation law is 


(8,2) 


— 


dx k dx l 

W d& kl ’ 


we call Atf the components of a co variant tensor of the 
second order. 

Further, if we have N 2 functions A\ whose transfor¬ 


mation law is 
( 8 . 8 ) 


A * ” dx* 13 


we call Aj the components of a mixed tensor of the 
second order. 

Note that the indices are placed on the tensors as super¬ 
scripts when they denote contra variance and as subscripts 
when they denote covariance. In particular, the mixed 
tensor A\ transforms like a contravariant vector with 
respect to the index i and like a covariant vector with 
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respect to the index /. Consequently the i is placed as a 
superscript whilst the j is placed as a subscript. 

Now choose A\ to be the Kronecker delta 6 k . From 
(8,3) we have 

j { dx l dzA h _ dx l dx k _ dx l _ 

5 dx k dx> 1 dx k dx* dx? * 


That is, the Kronecker delta is a mixed tensor of the 
second order whose components in any other system 
again form the Kronecker delta. This justifies the 
placing of one index as a subscript and the other as a 
superscript. Yet if we select the N 2 quantities d tf = 6) 
as the components of a co variant tensor in a coordinate 
system the components in the new system x* are 


- dx k dx k 
given by = 


Therefore the transformed com¬ 


ponents 6^ do not form the Kronecker delta. 


Ex. Prove that A^B^* is an invariant, if li i and C } are contra¬ 
variant vectors and a covariant tensor. 


§ 9. Higher order tensors 
A set of N*+* functions of the A 7 coordinates 

x i are said to be the components of a mixed tensor 
of the (s — p)-th order* contravariant of the s-th order 
and eovariant of the p-th order* if they transform ac¬ 
cording to the equation 


(811 ^ ^ Ah t_ . 

r i r s'■ - r u ‘ * ‘ dx** dx r i * 9x f j» 


Sic**! dx Q i 


, tn ■ * , 'f 1, 


on change of the coordinates x f to 0. This formula, 
although rather formidable in appearance, is merely a 
combination of (5.1) with respect to contravariant in¬ 
dices and of (6.1) with respect to co variant indices. 
The order of indices in a tensor is important. The 
tensor is not necessarily the same as the tensor A ji * 
(In the language of matrices A ji is the transpose of A a <) 
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If however two contravariant indices or two covariant 
indices can he interchanged without altering the tensor, 
it is said to be symmetric with respect to these indices. 
We shall now prove that if a tensor is symmetric with 
respect to two indices in any coordinate system, it 
remains symmetric with respect to these two indices in 
any other coordinate system* There is no loss in generality 
in proving this for the contravariant tensor — A**, 
On applying (8,1) we have 




as required* We cannot usually define symmetry with 
respect to two indices, of which one denotes contra- 
variance and the other covariance, because this symmetry 
may not be preserved after a coordinate transformation. 
The Kroneeker delta, however, is a mixed tensor which 
possesses symmetry with respect to its two indices. 

When all the indices of either a contravariant or a 
eovariant tensor can be interchanged without altering 
the tensor, it is said to be symmetric* A symmetric tensor 
of the second order has at most %N(N + 1} different com¬ 
ponents* 

A tensor each component of which alters in sign but not 
in magnitude when two contravariant indices or two co- 
variant indices arc interchanged is said to be skew- 
symmetric with respect to these indices. It can be 
shown by equations similar to (9*2) that the property 
of skew-symmetry is also independent of the choice of 
the coordinate system* Skew-symmetry, like symmetry, 
cannot be defined with respect to two indices, of which 
one denotes contra variance and the other covariance. 

If all the indices of a contravariant or a covariant 
tensor can be interchanged so that the tensor changes its 
sign at each interchange of a pair of indices, the tensor 
is said to be skew-symmetric. A skew-symmetric tensor 
A li of the second order has at most %N(N — I) different 
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arithmetical components, as all the quantities A u (no sum¬ 
mation) are zero* The several components of a skew- 
symmetric tensor of the A T -th order are either zero or differ 
merely in sign* So there is essentially only one non-zero 
component of such a tensor. 

The most important deduction from (9*1) is this: if all 
the components of a tensor in one coordinate system are 
zero at a point, they are all zero at this point in every 
coordinate system. Further, if the components arc iden¬ 
tically zero in one coordinate system, they are also iden¬ 
tically zero in every coordinate system* It is this property 
which constitutes the importance of tensors in physical 
applications. 

When a tensor is defined at all points of a curve or 
throughout the space V# itself, wc say that it constitutes 
a tensor-field. 


Ex. 1* If A is a skew-symmetric tensor, prove that 
+ S[^)A a = 0 . 

Ex, 2. Prove that the transformations of tensors form a group* 


§ 10. Addition, subtraction and multiplication of 
tensors 

It is clear that we cannot expect to give any tensorial 
meaning to the expression A i} + because it cannot 
satisfy the transformation law (9*1). It does, however, 
follow from this equation that any linear combination of 
tensors of the same type whose coefficients arc invariants 
is a tensor of the same type. For example, from the two 
tensors Aj k and Bj k , we can form the tensor XA) k + t*Bj k 
which will satisfy (9,1) provided that X and ft arc in¬ 
variants, in particular, A\ k + Bj k and Aj k — Bj h are 
called the sum and difference respectively of the two 
tensors. As another example, wc can write 

^a — i(^ij + ^fi) + — ^a)* 

Now A H 4- A h is symmetric and A iS — A H is skew-sym- 
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metric. Tims any covariant tensor of the second order 
is the sum of a symmetric and a skew-symmetric tensor. 
This is, of course, also true of a contravariant tensor of 
the second order. 

Let us select two tensors, one of contravariant order 
£ and covariant order p, the other of contravariant 
order i and covariant order q. It then follows from (9.1) 
that the products of the components form a mixed tensor 
of contravariant order $ + t and covariant order p + q t 
This tensor is called the outer product of the two 
tensors. For instance A^ nt — B** C l mnt is the outer pro* 
duct of the two tensors and C l mnt and is a tensor of 
the type indicated by its indices. 

The division, in the usual sense, of one tensor by another 
is not defined. 


§ 11. Contraction 

Let us start with any mixed tensor, say A^ mn , and form 
the sum A^ mr From (9.1) wo have 

j st dx^ d& dx l dx m dx n iS 

pqr " a? dx*dF» W* Almn ' 

Therefore 

j gr _ dx a dx r dx* 

VQT dx* Bx i dx p dx q dx T " lmn 
= cte™ 

dx i dx p dx q * lmn 
dx* d& dx m tn 
dx l dx p dx q “ lmn * 

Thus we observe that A is a mixed tensor, contra¬ 
variant of the first order and co variant of the second order. 
This process, which is called contraction, enables us to 
obtain a tensor of order r 2 from a mixed tensor of 
order r. In the above example we could contract a stage 
further and arrive at the covariant vector A When 
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contracting, any superscript may be used to sum with 
any subscript. Therefore we can form the following dif¬ 
ferent tensors by contraction:- A^ mp A \ 3 in , A*^ nf A^ mV 
At, A« mn , At A^, At, A*/ mi , At and A% n . If 
the tensor Af mn possesses any symmetric properties, 
there will be fewer tensors formed from it by con¬ 
traction. As another example, the invariant Aj is formed 
by contraction from the mixed tensor Aj. This justifies 
us in calling an invariant a tensor of zero order. 

We can also combine multiplication and contraction 
to produce new tensors. From the tensors A J? and 
we may obtain such tensors as A*J B^ nV B l inj , A£ B^ i{ 
and many others. This process is called inner multipli¬ 
cation of two tensors and the resulting tensor is called 
an inner product of the two tensors. 

Note carefully that we never contract two indices of 
the same type as the resulting sum is not necessarily a 
tensor. Also it should now be clear that with our index 
notation, the summation convention generally applies 
to two indices one of which is a superscript and the other 
a subscript, 

§ 12, Quotient law 

Sometimes it is necessary to ascertain whether a set 
of functions form the components of a tensor. The direct 
method requires us to find out if they satisfy a tensor 
transformation equation of the type (9.1). Li practice 
this is troublesome and a simpler test is provided by the 
quotient law. The quotient law states that N p functions 
of x i form the components of a tensor of order p s ( whose 
contravariant and covariant character can readily be 
determined), provided that an inner product of these 
functions with an arbitrary tensor is itself a tensor. It 
will suffice to set out the proof for the following parti¬ 
cular case. The set of N 3 functions A m form the com¬ 
ponents of a tensor of the type indicated by its indices if 
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= C*\ 

provided that Bf f is an arbitrary tensor and C pk a tensor* 
The transformed quantities, referred to a system of coor¬ 
dinates x\ satisfy the equations 

A^Bl = 

which become on substitution from (0,1) 


Sm dj? factor n t 

dx l dx 1 30 mn 


d& 30 
dx q dx f 


dx p dx k 

zzi _ 4 Ur m 
dxf 1 dx f H * 


With a change of dummy indices, we have 


dx 7 * 

dx 1 


\ A m d jf_ _ Am «r ?i*l 

L dx* dx j dxf J 


Kn 


= 0. 


On multiplying this equation by dx*/dx p and summing 
over p from 1 to N, (in future we shall merely write *on 
inner multiplication by dxfjdx p ') we obtain 


( 12 * 1 ) 



d : v^dx»_ Afnnf 05*1 
d& d& 3;c r J 


= o. 


Since is an arbitrary tensor, we can arrange that 
only one of its components differs from zero* Now each 
component of B^ n may be selected in turn as that one 
which does not vanish* This shows that the expression 
in brackets is identically zero. That is 


( 12 * 2 ) 


Am dx™ dor _ 30 

dx* 30 dx T ' 


■ * * » * dxf dx* 

Inner multiplication of this equation by yields 

the result 


dx* dx 1 dx k 
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Thus A mnr is a tensor of the third order and contra- 
variant in all its indices* In the above proof it is important 
that the tensor shall be arbitrary' and must not 
possess any symmetric or skew-symmetric properties* Let 
us examine what happens if is symmetrical in m 
and 7L It is no longer possible to deduce (12*2) from (12*1). 
The correct deduction now is that 


Si » ^ , A*’* — — - A™ — 4 A™ — 

dx* 30 30 30 “ 30 1 30 ' 


On changing several dummy indices, this equation be¬ 
comes 


+ &*) 3 ^ 30 = ( Amnr + Anmr ) 


dx h 

dxf ' 


Inner multiplication by ^ shows that A mnr + A nmr 

is a contra variant tensor of the third order. If in 
addition we know that A mnr is symmetric in m and n, 
it follows immediately that A mnT is a symmetric tensor 
with respect to the indices m and n . We observe, there¬ 
fore, that the quotient law must be applied with care* 

Ex* 1* If ^1' and 7J* are arbitrary contra variant vectors and 
C ii A i B i is an invariant, show that C ff is a covariant tensor 
of the second order. 

Ex, 2. If A * is an arbitrary contra variant vector and is 

an invariant, show that C iS -p C lt is a covariant tensor of 
the second order* 


§13* Conjugate symmetric tensors of the second 
order 

Consider a symmetric covariant tensor of the second 
order A a whose determinant \A U | ^ 0. Let B'* denote 
the expression formed by dividing the cofaetor of Ay in 
the determinant | A if | by \A (i | itself* We shall prove 
that the B ij so obtained are the components of a eontra- 
v &riant tensor of the second order. In anticipation of 
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this, B 1 * is labelled as if it were a contravariant tensor* 
We have from the theory of determinants** 

(13*1) Atf B ik = <5J. 

We cannot establish the tensor character of B ik by 
applying the quotient law directly to this equation, 
because A u is not arbitrary* Let us choose an arbitrary 
contravariant vector CK Then D t = A ti O is an arbitrary 
covariant vector, because these Nequ a Hons can be unique¬ 
ly solved for the 0 in terms of the since \ A iS \ ^0. 
Consequently, 

D f B ik = A u a B** = = CK 

Now if we apply the quotient law to the equation 
D f B ik = C k , we see that B ik is a contravariant tensor of 
the second order. Also it is clear from the definition that 
B 1 *, like A w is symmetric. 

We will next attempt by the same process to obtain 
another tensor from B ik . Let E denote the cofactor 
of B ij in the determinant | B i} | divided by | B if | itself. 
From the theory of determinants | A if | * | B u | = 1 and 
consequently | B i} | ^ 0, which means that E ti always 
exists. Further we have 

E is B ik = 3J. 

Inner multiplication by A ki yields on application of 
(13,1) that 

Eij = A n = A u * 

Thus this process only leads back to the original covariant 
tensor of the second order* We say that A u and B li are 
conjugate tensors if they satisfy equations (13.1), It is 
important to note that a tensor of the second order has 
a conjugate only if its determinant is not zero* 

Ex, If A^ = 0 for i show that the conjugate tensor B** ^ 0 
for i ^j, and B u — 1 jA u (no summation). 


* A, C, Ait ken, Determinants and Matrices, pp* 51—52* 
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§ 14, Fundamental tensor 

At this stage we introduce the concept of distance into 
our space V If the distance ds between the neigh¬ 
bouring points with coordinates x* and x* + dx t is given 
by the quadratic differential form 

(14.1) d$ 2 = g ii dx i dx i 

where the g if are functions of x f , subject only to the 
restriction g — | g 0 | ^ 0, the space is said to be a 
Riemannian space* In addition we postulate that the 
distance between two neighbouring points is independent 
of the coordinate system. That is, ds is an invariant* 
From the quotient law, it follows that g^ + g H is a co¬ 
variant tensor of the second order. We can write 

in ~ + Bn) + i(Sif — Bn)- 

The contribution of |(g f , — g ii )dw i dx i to ds 2 is zero, hence 
there is no loss of generality in assuming that g ti is sym¬ 
metric, Tims g ti is a covariant symmetrical tensor of the 
second order called the fundamental tensor of the Rie- 
mannian space* The quadratic form g ii dx i dx i is called the 
metric* It is also the square of the line-element ds. 
The line-element ds of a three dimensional Euclidean 
space, referred to a system of rectangular cartesian axes, 
is 

ds 2 = (da; 1 ) 2 + (tfc 2 ) 2 + (cfo 3 ) 3 * 

All the components of the fundamental tensor are zero 
except g n = g 33 = g 33 = i. It is evident that the metric 

C 




18 


TENSOR CALCULUS 


§ 15 

of a Euclidean space is positive-definite*. That is, da 2 
is zero when d# 1 = dx % = dx* = 0 but can only take 
positive values for all other real values of da? 1 , dx 2 and cfe a . 
The Special Theory of Relativity discusses the four 
dimensional space with the line-element ds given by 

(14.2) ds® = — (dar 1 ) 2 — (da? 2 ) 2 — (tfe 3 ) 2 + 

This metric is not positive-definite, as it is positive for 
all curves along which a 4 , x % and # 3 are all constants, but 
negative for all curves along which is constant. Thus 
along these latter curves the distance between neigh¬ 
bouring points cannot be real In order that the distance 
ds between two neighbouring points be real, equation 

(14.1) will be amended to 

(14.3) ds a = eg\ ij dx i dx i t 

where the factor e called the indicator takes the value 
+ 1 or — 1 so that ds 2 is always positive. 

Ex* Show that the metric of a Euclidean space, referred to spherical 
polar coordinates x 1 = r, x* = 6 and = ip is given by 
= dr 1 -f r 2 d0 2 + r 2 sin a Ody** 

§ 15. Length of a curve 

Consider the curve w* — x*(t) with the parameter t. 
From (14.3) the length of the curve between the points 
corresponding to t = t t and i = f s is given by 

V % \I d® 1 dot j. 

(15.1) '-J,, I «*![*■*• 

dx^ dx^ 

If g tj ^ ~ = 0 along a curve, then the two points 

corresponding to ^ and t % are at zero distance from one 
another, despite the fact that they are not coincident. 

* A. J* Me. Connell, Absolute Differential Calculus, p. 1G. 

A. C. Aitken, Determinants and Matrices, p. 137. 
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Such a curve is called minimal or null. The curve 
given by 

( x l = c \ r cos 6 cos ip dt> x* = c \ r cos 0 sin ip dt , 

05 * 2 ) J J 

[ x* = c j r sin 0 dt, ar 4 « J rdf, 

where r, 0 and ip are functions of t is a real null curve 
in the V 4 whose metric is (14.2), It is clear that no real 
null curves lie in a Ricmannian space whose metric is 
positive-definite. 

A curve will consist of portions along which the in¬ 
dicator e is + 1, portions along which the indicator is 
— 1, and null portions. The length of the curve is then 
the sum of the lengths of these portions, the null part 
contributing zero to the value of the length. 

Except in the case of null curves, the parameter t may 
be chosen as the are-distance $ from some fixed point 
of the curve. From (14.3) 

dx 1 dx s 

(1 »-3) gi, - « 

along any portion of a curve which is not null. 


§ 16. Magnitude of a vector 
The magnitude A of the contra variant vector A i is 
defined by 

(16.1) {Af = e w gttA*AU 

where e^ is the indicator + I or — 1 which makes A 
real The magnitude A is an invariant. In a Euclidean 
space, referred to rectangular cartesian coordinates, 

(16.1) reduces to the familiar definition of the magnitude 
of a vector. 

At this stage it is necessary to introduce the contra- 
variant tensor conjugate to g w which can conveniently 
be written g y . Then equations (13.1) read 

(16.2) g ijg * = <5*. 
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We can now define the magnitude B of the covariant 
vector B ( by the equation 

(1(5,3) 

where e {B ) is the indicator of the vector B if and it is clear 
that B is an invariant, 

A vector, whose magnitude is unity, is called a unit 
vector. It follows from (15.8) that dx { lds is a unit con- 
travariant vector* If the magnitude of a vector is zero, 
it is called a null vector* The tangent vector to a nut! 
curve is a null vector. 


§ 17. Associate tensors 

The inner product of the fundamental tensor g ti and 
the contravariant vector A i is the covariant vector g^A*, 
which is said to be associate to AK We define 

(17.1) A t = gii AK 

Similarly we define 

B' = gUB; 

and say that the vector B* is associate to the vector B r 
The relation between a vector and its associate is 
reciprocal, for the vector associate to Ai is 

g^A, = pg^A* = &*A* = AK 

This process of association is often referred to as "lowering 
the superscript’ or "raising the subscript 5 respectively. 
We have 

e ld) g"AW = e w g ii g ik A k g‘ l A l = e w g a A k A» 

which shows that the magnitudes of associate vectors 
are equal* 

Tfic process of raising and lowering indices can be perfor¬ 
med on tensors* From the tensor A^* m wc can form associ- 
ate tensors like A r ik tm =g ri A^ k tm or A i f’ t =g ir g u g mt A iit lm . 
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The dot notation is introduced to indicate which indices 
have been raised or lowered. The dots will be omitted 
when there is no possibility of confusion. For example, we 
shall write A f * = g iT g Ss A r9 . Note very carefully that 
although gii and g ti arc conjugate tensors, the tensors 
A if and A** are not as a rule conjugate. 


Ex, Show that {A) % = €(A\A t AK 


§ 18. Angle between two vectors - orthogonality 
The angle between two unit vectors A 1 and B l is defined 
by 

(18,1) cos 0 = g if A*B* = A S B* = g*A s B k - A*B k . 

This equation will be found to agree with the usual for¬ 
mula cos 0 =a IV + mu' + nri for the angle between the 
unit vectors (J, m, n) and (l\ m\ n r ) in a three dimensional 
Euclidean space when it is referred to rectangular car¬ 
tesian coordinates. We will now show that (18,1) always 
defines a real angle between two real vectors if the metric 
of the Riemannian space is positive-definite. For then the 
magnitude of the vector XA* + pB* is greater than or 
equal to zero for all real values of X and /*, That is 

g "(XA* + fiB i ){XA i + pB*) ^ 0, 
which reduces to 

A* + 2 Xft cos 0 + ft 2 S 0. 

That is 

(A + fi cos 0)2 + fi 2 (l — cos* 0) ^ 0* 

Since this equation is true for all values of A and pt, 
it follows that 1 —cos 2 0^0, Thus | cos 0 | % 1 which 
means that 0 is real. If the metric is not positive-definite 
then the angle between two real unit vectors need not 
be real. 

An immediate deduction from (18,1) is that the angle 0 
between two vectors A* and BK which are not necessarily 
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unit vectors, is given by 


(18.2) 


cos 6 = 


g hAW 

Ve [ A ) e {B) g lm A l A m g rs B r B>' 


Two vectors arc said to be orthogonal to one another 
if the angle between them is a right angle. From (18.2) 
the necessary and sufficient condition for the ortho¬ 
gonality of the two vectors A 1 and B* is that 
( 18 . 3 ) 0 . 

We do not define the angle between two vectors if one 
or both of them happens to be a null vector, but we shall 
take (18.3) as the definition of orthogonality of two null 
vectors. It follows that a null vector is self-orthogonal. 


Ex. Prove that (1, 0, 0, 0) and (y% 0> 0, ^B/c) are unit vectors in 
the f 7 4 with the metric (14.2). Show also that the angle between 
these vectors is not real. 


| 19, Principal directions 

From the symmetric covariant tensor A ti we can 
construct the determinantal equation 

(19*1) | Atf — Xgn | = 0, 

which is of degree N in A, When we change to a new 
coordinate system x\ this equation transforms to 


— ^gtk) 


dx 1 dx k 
Bx i dx* 


= 0. 


On applying the multiplication law of determinants, this 
equation can be written 


A ik — hg lk 


dx 1 
dx i 


® 0 . 


The Jacobian determinant 


dx 1 


dx* 


sequently this equation reduces 
I AIk — %gik \ 


does not vanish; eon- 
to 

= 0, 
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By comparison with (19*1), we deduce that the roots 
Xtf) of this equation are invariants. (Capital letters are 
used to label the roots. Their enclosure in brackets em¬ 
phasises that they have no tensorial significance, and the 
summation convention is not to apply to them.) 

Now consider the N equations 

(10,2) gff)£(fQ — A* 

where X [K) is a simple root of (19.1). These determine 
the ratios of the N values of L\ K) . We cannot deduce 
immediately from the quotient law that L\k) is a contra- 
variant vector, because the tensor — X(?q ga, with 
which it shares inner multiplication, is not arbitrary. 
Instead we change to the coordinate system x*, and equa¬ 
tions (19.2) transform into 

(A, k — Wit) = °' 


Inner multiplication by dx*jdx m yields 

dx 1 

(Ai m — hlQilm) L \tq = 


These N equations determine the ratios of the N quan¬ 
go-i 

tities ^which arc the components of in the 

0 coordinate system. That is, h \ K j transforms in accor¬ 
dance with equations (5.1) and thus Lj^j is a contra- 
variant vector. Let us now choose the ratios of the Lj jq 
so that it is a unit vector. That is, 

(19,8) q — 

where e^ is the indicator of the vector L|jq. 

Similarly we can show that corresponding to each 
simple root A(jf) of equation (19.1) there corresponds the 
unit vector satisfying the equations 
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and 

^ 15 ^ J^( M) &{Af j ) — £{Af) - 

When A ( mj is a multiple root of (19.1), the vector L‘ 
is not uniquely determined by equations (19.4) and (19.5). 

Let us choose two simple roots X {K) and X m of (19.1), 
Smee X (K) ^ A (W) , the inner multiplication of (19.2) by 

L\ U) and (19.4) by L\ K) gives us on subtraction, the 
equation 

This shows that the two unit vectors L* and L\ m are 
orthogonal. Thus if all the roots of equation (19.1) are 
simple at any point, the covariant symmetric tensor deter¬ 
mines uniquely N mutually orthogonal unit vectors. The 
directions of these vectors at a point are called the 
principal directions determined by A i} . If the metric 
g M dvW is positive-definite all the roots of (19.1) are 
real*. That is, the principal directions are real in a space 
with positive-definite metric. 

Now consider the invariant X defined by 


3 _ A„L*V 

The finite maxima and minima of X are (riven bv 
that is, by y dL 1 

A u^i m L l V») — gii L*(A lm L l V”) = 0 . 
This equation can be written 


= 0, 


— = o. 

Elimination of the TJ yields the determinantal equation 

Mo — Aft, I = 0. 

Thus the finite maxima and minima of X are those values 
corresponding to the principal directions determined 
by A u- 


* W. L. Ferrar, Algebra, p. 145. 
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If Afj = Xg u at a point, then the principal directions 
are indeterminate at that point. If A (j = Xg i} at all points 
of a V N , the space is said to be homogeneous with 
respect to the tensor A i} . 

In a Euclidean space of N dimensions referred to 
rectangular cartesian coordinates, the components of the 
fundamental tensor g if form the unit matrix. Hence the 
roots of (19.1) are in this case the latent roots* of the 
matrix A (i and the principal directions are those of the 
latent vectors. 

Ex. Prove that A (t ^‘\ k) ~ e (K) ^-(X) an< ^ A v ^\k) B (Jf) “ 
(No summation over {£)). 

Examples 

1. Show that the angle 0 between the vectors A* and B‘ is given by 

Sin s g _ ( e (/i )*(£)&<&*—71* . 

«(B) ghigikA l A i IPB t 

2. If An is a skew-symmetric covariant tensor, prove that (/I,,/ Vg, 
A (I /V£ AjVg) arc the components of a contra variant 
vector. Show also that (Vg.4* 5 , v'g/1 31 , VgA li ) are the com- 
ponents of a co variant vector if A** is a skew-symmetric 
contravariant tensor. 

3* Prove tliat no relation of the type 

tytistkk + -f v£ ik A hi — 0 

can exist in a Ftf, {N >1), where A, /* and v arc invariants 
and is a symmetric tensor. 

If A hk is a skew-symmetric tensor, show that this equation can 
not exist in a for which N > 2. 


* A. C. Aitken, Determinants and matrices, p. 78. 
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COVARIANT DIFFERENTIATION 

§ 20* Christoffel symbols 

Although we found in section 5 that dx i jdn is always a 
eontravariant vector, the exercise at the end of that sec¬ 
tion should convince us that its derivatives dW/du® do not 
form a vector, whose components in any other system are 
the corresponding second derivatives* Again, in section 6 
we proved that the partial derivatives of an invariant 
form the components of a eovariant vector; yet we shall 
show in section 22 that the derivatives of a vector do 
not form a tensor whose components in any other system 
are the corresponding derivatives of the transformed 
vector* Our aim is now to build up expressions involving 
the derivatives of a tensor, which arc the components of 
a tensor* In order to carry out this programme we must 
first investigate two functions formed from the fun¬ 
damental tensor g ijm These are the Christoffel symbols of 
the first and second kinds defined respectively by 

<*«■•> + 

and 

( 20 . 2 ) 

Although we shall see that the symbols [£j, k] and 

are not tensors, this notation is introduced in keeping 
with the summation convention which generally applies 
to two indices, one a superscript and the other a subscript* 
Thus all but one of the indices of the Christoffel symbols 


i) 
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are regarded as subscripts. The exception is the l in the 
symbol of the second kind, which is treated as a super¬ 
script* 

The definitions show that both symbols are symmetric 
with respect to the indices i and j. Inner multiplication 
of {20.2) by g tm yields 

(20.8) [*?, m] = g, m | t ?.| • 

It follows immediately from (20.1) that 

(20.4) ^ = [ij, k] + [kj, q. 


We now wish to express the derivatives of g** in terms 
of the Christoffel symbols, and so we differentiate equa¬ 
tion {16.2) with respect to a; 1 , and obtain 


dx l g,i h dx' 8 


Inner multiplication by g im gives us 

-L g Srngik — 0. 

dx 1 ^ 8 8 dx 1 

Substituting from (20.4) and (20*2), we finally obtain 

< 20 ' 5 > W -«“ («)-«“{”)• 

We shall now deduce a useful expression for j^j* Diffe¬ 
rentiate the determinant g = \ g tj remembering that 
g* m g is the cofaetor of in this determinant, and obtain 

dtf * * fat 

From (20,1) and (20*2) and the symmetry of we have 









28 


TENSOR CALCULUS 


§ 21 


f M = I gin, f^km + d lm_d2tj\ J_ im 

\ i ?i 2 s a®* a^ m j 2 s g**' 


Thus 

( 20 . 0 ) 


f n i i ag a . /- 

Uj = 2 73S = ai5 (lo ® v * > - 


Since g is not an invariant, it does not follow that 
| „} is a eovariant vector. If g is negative, equation (20.6) 
should be altered to {.*.} = A {log V^}. 


Ex. 1. Calculate tlie Christoffel symbols corresponding to the 
metrics 

(a) = (<&*)* + (s 1 )*^ 8 ) 8 + {t*rsln*x*(&a)K 

(b) d$* = (dri) 1 -j- (?(£*, a^Hdk 3 ) 8 , where G is a function of 
x l and x\ 


Ex. 2* If the metric of a V N is such that g iS = 0 for i ^ show 
that 

f i 


.0, (M—_ 

U2 J Hu 


where i s j and k are not equal, and the summation conven¬ 
tion does not apply. 


§ 21. Transformation law of Christoffel symbols 

The fundamental tensor g t) , being covariant, trans- 
forms according to the equation 


( 21 . 1 ) 


, _ dx* dx* 

gtm ~ d&d&» g *** 


On differentiating this equation with respect to x n \ we 
have 

dgim fa* dgti dx* ( d\x* dx* dx i d 2 x* 

dz n dx l dx m dx k dz n 1 dx l dx n dUF* ^ d£ l 
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We subtract this equation from the sum of the two 
similar equations obtained by cyclic interchange of the 
indices l, m and n, and divide by two. Then with appro¬ 
priate changes of dummy indices we have 


r!-, ,, 3** dx* dx k , _ dx* d a x i 

(21.2) [Im, n] — [tj, k ] ^ d£m g. n + go g £fl * 


where the bar over the Christoffel symbol indicates that 
it is calculated in the coordinate system x l with respect 
to its fundamental tensor g fJ * Now the transformation law 
of the contravariant fundamental tensor is 

dx n dx p 

(si.8) 


Inner multiplication of both sides of equations (21.2) by 
the corresponding sides of (21.3) yields on reduction the 
relation 

j~p 1 (s \ dx p dx l das* , dx? d 2 £* 

[(m j daf dx m ^ dx J dx l dx^ 

Equations (21*2) and (21.4) comprise the transfor¬ 
mation laws of the Christoffel symbols, and clearly in¬ 
dicate that they are not tensors* However, in the very 
special case of linear transformations of coordinates, 
d^&tjdxldx™ = 0 and the symbols then transform like ten¬ 
sors. Now inner multiplication of (21*4) by dufjdx* yields 

9 2 ^ f _ f V 1 f r\d& dx^ 

dx l dx m — \imj Bx p dx l dx m 

This important equation expresses the second partial 
derivatives of the of with respect to the z* in terms of 
the first derivatives and Christoffel symbols of the 
second kind. 


Ex. Prove that the transformations of Christoffel symbols form 
a group. 
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§ 22. Covariant differentiation of vectors 
Let us investigate the tensor character, if any, of the 
partial derivatives of a contra variant vector. We start by 
differentiating the transformation law 


( 22 * 1 ) 


A* — A‘ — 


with respect to x and obtain 

dA k _ dA‘ dx" dx k t d 2 x k dx" 
dx* dx n dx* dx 1 1 dS;* dx n dx* 


The presence of the last term on the right-hand side of 
this equation shows that the partial derivatives dA k jdx> do 
not form a tensor. To obtain a tensor involving the 
partial derivatives we eliminate the partial derivatives 
of the second order by means of (21.5) and this gives us 

dA k = dA* d& ete* - f|yi dx*_fk\dx T a*»"| 

dx 1 dx" dx i dx* " r J dx 1 |_(in j 5x» j d£‘ d£"J 

In virtue of (22.1) and by appropriate changes of dummy 
indices, this equation reduces to 

M* _L f *1 AT __ [ dA ‘ , i T I *rl dx" dx k 

dx 1 ( r/j \_dx n ~ r (raj J dx 1 ' 

Introduce the comma notation 



so that the above equation can be written 

4>c _ Xi dx k 
- 1 ’ a dx‘dx‘* 

Henec from (8.8) it is obvious that A k } is a mixed tensor 
of the second order and it is called the covariant deri- 
vative of A k with respect to xK 
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In order to set up the corresponding entity for covariant 
vectors, we may conveniently start by differentiating the 
transformation law 


(22.3) 



with respect to xK This yields 

3 At _ 3 Aj 3 ^ , , 3 2 ^ 

dx 1 — dx k dx 1 dx* i dx l dx l 


Again, by (21.5), we eliminate the partial derivatives of 
the second order. Further we change dummy indices as 
required and we have on substitution from (22.3), that 

dtifrXx == \dJ 1 _ir) A ld^d^ i 
dx 1 \ il) A m [<*•" b'n| f J dx* dx 1 

The comma notation 



is now introduced, and the above equation can be written 
j dx* 3# n 

Ai > t = A, ' n W i dJ’ 

which shows that A Un is a covariant tensor of the second 
order, called the co variant derivative of A, with res- 
peet to x n * 

In a Euclidean space of N dimensions, referred to 
rectangular cartesian coordinates, the components of the 
fundamental tensor are zero except g n — g 22 = . . * — 
— g NN = i, Thus all the Christoff el symbols are zero and 
so covariant differentiation reduces to the familiar partial 
differentiation* It is well to observe that the Chnstoffel 
symbols do not all vanish in a Euclidean space referred, 
for example, to spherical polar coordinates* 

We can construct the invariant A*f by contraction* 
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Applying (20.6) we obtain 


A - i = W + [lj\ A ' = W + A * 0°& ^sb 


That is 
{22.5) 


Ai -* sTg da? ^ At )‘ 


This invariant is called the divergence of the contra- 
variant vector A i and is often denoted by div A*, Tiic 
divergence of a covariant vector A t is defined by 

(22,6) dr yA i =#*A Stk . 

Tiie partial derivatives of an invariant form the com¬ 
ponents of a covariant vector. We extend the definition 
of covariant differentiation to invariants by calling the 
familiar partial derivative the covariant derivative. That 
is, from the invariant I we form by definition 


X.i 



Since I t is a covariant vector, we find from (22 A) that 
its covariant derivative with respect to x i is given by 


(22,7) 


(*«>,* 


d*I 

dx*dx* | ijf dx r * 


Thus (/ *)/ — (1 ^ That is, the covariant differen¬ 

tiation of invariants' is commutative. In section 31, we 
shall see that the covariant differentiation of vectors is 
not in general commutative. 

We can form the divergence of the covariant vector 
I this is called the Laplacian of / and written V 2 L Then 


**(*.*),* = 


i 91 l 

’M£\ 

\dx*dx k \ 

jkj daf) 


Ex. 1. Show that — A ntj = 

s*** 


ix> 
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Ex. 2. 


Ex. 8, 


Prove that div A f = 


V i u X* 


Evaluate div A* and V % 1 in (a) cylindrical polar, (b) spherical 
polar coordinates. 


§ 23, Covariant differentiation of tensors 

In the last section we constructed tensors containing 
the partial derivatives of vectors. Can we now extend the 
process of covariant differentiation to tensors? For our 
typical tensor we shall choose There will be no loss 
in generality involved as this tensor has one contra- 
variant and one covariant index. Inner multiplication of 
its transformation law (8.3) by dx m ld& yields 


(23.1) 


Ai djr 

Ai dx‘ 


_ dm 

— 


dx l 
dx s * 


We differentiate with respect to x k 9 then eliminate the 
second order partial derivatives by means of (21,5) and 
obtain 


dAj dx m , j t H p 1 dx m (m\ daf dx*l 

dx k dx* 1 * dx p \rsj 9.x 1 dx k j 

_ dA? 4m r/ P1 — f M — 

W W k d& + 1 Ll**j l«j d& dx k j ‘ 




Applying (23,1) to this and clianging appropriate dummy 
indices, this equation becomes 

rMI i A n l T )-A f R1 — 

Las* U) MwJ a** 

~ldx t ' 1 \r*j r IfljJ dx* d& 

Introduce the comma notation 

(“■ 2 > + 


D 
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and the above equation takes the form 
Ti 8x™ _ dx* dx l 

*' k T& ~ Al % di* W' 

On inner multiplication by dx T ldx m we have 
Jr _ 4m d ® 1 

i r h ^ g-fc * 

Hence Af t is a tensor of the third order of the type 
indicated by its indices. This tensor is called the eo^ 
variant derivative of Af with respect to 
An examination of equation (28.2) shows that the 
eovariant derivative of Af contains three terms. They 
are (1) the partial derivative, (2) a term with positive 
sign similar to that contained in the covariant derivative 
of a contravananf vector and (3) a term with negative 
sign similar to that contained in the covariant derivative 
of a covariant vector. This suggests that the expression 


(28.8} 




,u t 

r P* n 



+ T 

d-i 




is a tensor, called the eovariant derivative of A™* 
with respect to of 1 * The proof which we worked out* for 
the mixed tensor of the second order applies also to (23,8). 
We will omit its tedious details and supply an alternative 
proof in section 80. 

On referring to equations (20,3), (20.4), (20.5) and 
(23.3) we deduce that 


(23,4) 



(23.5) 
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S 24 
and 
(23.6) 




The eovariant derivatives of eovariant derivatives arc 
again tensors, Wc indicate these eovariant derivatives 
of the second order by adding another index without a 
comma. For example A Uk is the covariant derivative of 
A iif with respect to x k . 

Ex- 1. K A if = B ui — B f , „ prove that k 4- A fk , ( + A kUj — 0. 
Ex* 2. By means of (23,5), show that div A* = div A t , 

Ex, 8, If A iik is a skew-symmetric tensor, show that 

——— {^g A 4 **) is a tensor, 

VB ^ 


§ 24. Laws of covariant differentiation 

Covariant derivatives obey the following laws:- 

(1) the covariant derivative of the sum (or difference) 
of two tensors is the sum (or difference) of their covariant 
derivatives. This law is an immediate deduction from 
(23.3), 

(2) the eovariant derivative of an outer (or inner) 
product of two tensors is equal to the sum of the two 
terms obtained by outer (or inner) multiplication of each 
tensor with the covariant derivative of the other tensor. 
Consider as an illustration 

-^ A ^+ A « (A )•>' [{AH*+(AH 

= A ihm B' + A ti B\ m . 

This type of proof is quite general and will apply to any 
ease of outer multiplication of two tensors, (Another 
proof is provided in section 28). Contraction of l and j 
gives us 

{A il B*U = A ii . m B’ + A it & m , 
which makes it clear that the rule also applies to inner 
products. 
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{3) the tensors g w , g* and dj are constants with respect 
to covariant differentiation. This is merely another way 
of stating equations (23,4), (23.5) and (23,0)* For example 

{^An\ m = g^An + g**A € Um — g { *A fltm . 

Ex* If 1 and J are invariants, show that 
div = 


§ 25* Intrinsic derivatives 
Consider the tensor A^y^* whose components are 
functions of t along a curve V = x l {t). The intrinsic 
derivative is defined by 


(25*1) 


Kf::5 

to 


= 


r„lc 


dx k 

It 


Accordingly, the intrinsic derivative is a tensor of the 
same order and type as the original tensor. 
Corresponding to the invariant I, we have 

61 dx* _ dl_<kt _dl_ 
dt '* dt dx k dt ~ dt 


That is, the intrinsic derivative of an invariant coincides 
with its total derivative. 

Intrinsic derivatives of higher order are easily defined. 
For example. 


_ 8 (6Afi / t dx l 

~8fi' - to \li ~) ~ \ i‘* It ) tl ~dt ' 


In general, intrinsic differentiation is not commutative. 

From (22.2), (22.4), (23.4), (28.5) and (23.0) we cal¬ 
culate that 


(25.2) 

(25.3) 


dA* _ dA* , 

fjfcl 

1 A 

dt ~ dt _r 1 

WJ 

\ A, w 

dA k _ dA k I 

frl 

1 A ^ 

dt ~ dt \ 

tf\ 

| A *li 
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and 

(25.4) 


fa, = W = W 

dt dt 6t 


From the definition of intrinsic derivatives, it follows 
that they obey the same three laws which apply to 
covariant derivatives. 


Ex, Show that 


S /daA _ ( J i \ da* dri* 

<5f \ dt) = ~dF + \ j*J ~di ~di 


p, 2S Ex. 1 
kind arc 

(a) 


Solutions 

Tlie only non-zero Chris toff el symbols of the second 


l/* l t 


-" l> (as)- (m) 

= — gin ** cos **, ( * 1 = l/a l , 1*1™ cot*’; 

_ ia<? f 2 l = 1 2 1 = -L^ 

|l2j^2Gba:'’ |22j 2Cbs s ‘ 


(b) 

p* 33 Ex* 3* 

, . .. .. b.4 1 , M* 2U4* 1 

(a) dlVj4 + ^ + 

+ _LMi + !d: + l Ati 

lx 1 tto 2 bar 3 X 1 


wn -= 


b*/ 




b*2 


bs* 

(b) div A* — -f —- + —- + — A 1 + cot x*A*, 
1 ba: 9 ae 1 


0 #)* (x'Ffb&y 0®*) 1 


bar 1 

%A t 


1 *A t 

ba: 1 + (zM* bs* + 


cotg* i . 

(a 1 ) 1 sin a x 1 ba: ! ^ a: 1 1 (x ! ) 4 




b»i 


b 1 / 


2 b/ cot x* bJ 


"(bar 1 ) 1 (x 1 )* (bar 3 )* ' (a^sin^tb,**)* a: 1 bar 1 (a: 1 ) 1 bar* 

The corresponding results in Rutherford's Vector Methods, pp* 
72-^-73 t differ slightly because they arc expressed in terms of the 
physical components of the above vector (see section 62)* 


























CHAPTER IV 


GEODESICS - PARALLELISM 
§ 26* Geodesics 

In Euclidean three-dimensional space, a straight line 
is the path of shortest distance between two points, and 
it is our aim to generalise this fundamental concept to 
Riemannian spaces. Let C be the curve x* = x*(t) with 
parameter t joining two fixed points P 0 and P v whose 
parameters are * 0 and t x respectively. Then the distance 
$ along the curve between P 0 and P x is given by 

(2 “> -£)£??* 

Consider all the curves passing through the two fixed 
points Pq and P v Any of these curves, for which the 
distance f J 0 P 1 measured along the curve is stationary, 
is called a geodesic. We could obtain the differential 
equations of the geodesics by applying Euler’s equations, 
a well-known result in the calculus of variations, to 
equations (2C.1). However we shall find it more instruc¬ 
tive to appeal to first principles. 

Let us choose a small arbitrary vector dx* varying con¬ 
tinuously along C. Then the equations x* = x* + dx* 
define a neighbouring curve C to C. Further let us impose 
the conditions that = 0 at P 0 and P r This means 
that the curve C always joins P 0 to P v The distance s 
from P Q to Pj along C is given by 
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where the g H (£) are now functions of x l . We have 

.. d{dx*)\ 


._.dx*dx* ( dg„ \ (dx* d(<5^)\ (dx 

8 ‘» IT rfT=l 8 " + a ?") v w + sr) u 


dx*dx* 

= g(i ~dt~^ +gii 


/dx* 



1 dt , 

dx*d(dx*) , i 

dt 

dt 


1-3T+- 


at } 

Hu . „ dx* dx* 
a®* 5 dt dt ’ 


where terms of higher order than the first are neglected. 
Thus 


Y 

Y 


dx* dx* 

e ^ x) ~ZTi- 


' 8 »rfT 


dx* dx* 
dt 


1 + 


dx* d(6x*) , 1 dg if . k dx* dx*~l 

dx* dxi 
gii dt dt 


Consequently the variation in length ds from the curve 
C to the curve C is given by 


= s — s 

rh 


dx 4 d(dx*) 1 dg iS A 

^ Hi dt 2 dx k dt dt 


i/ dx* dx 4 

r egii ~di!u 


dt. 


We now simplify this equation by choosing the arc- 
distance s along C as parameter and obtain 


3s 


f*i f dx* 

-J, 


d{dx*) , 1 3g it dri dx? 

ds i 2 3ai* ds ds 


ds. 


where s 0 and s t are the values of s corresponding to the 
points P 0 and P x respectively. Integration by parts yields 

. f dxi t,, r,. v d ( dx<\ 13 g ik dx'dx*-] 

*• - L 8 " * & Ir J.. *r* l 8u nr* ® s sr 

The integrated portion vanishes since 6x* is zero at P 0 
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and P v Also we have 

±( <fa* 

ds \ gi * ds)~ gi * ds* + 3^"* 


§ 26 


1^ dg^ Ar* dx k , I 9g w da:* 
^ ds 2 2 3a? 12 <£s ds 2 3a:* ds ds 


It follows that 


for o\ ju f* 1 * ,r . r*f ^dx l dx k 1 

( 20 . 2 ) « s —a*< [g„ w + [.*,,] *. 

The variations 8x f are arbitrary, thus the necessary and 
sufficient conditions that the curve C be a geodesic arc 

(26.8) g W ^r+ ^,7]^-^-=°. 


Inner multiplication by g Jl yields the contravariant form 

(26.4) - (—) = ^£.‘ 4 - I M — — - 0 

1 - ds 2 + (rt) ds ds - °- 

Either set of equations (26.3) or (26.4) are the differen¬ 
tial equations of a geodesic. They constitute N differen¬ 
tial equations of the second order. The theory of differen¬ 
tial equations states that a solution — *'(s) is deter¬ 
mined uniquely if the initial values of x f and dx i /d$ are 
given at any point. Geometrically this means that there 
is a unique geodesic with given direction at any point 
of the space. We defined the geodesic in terms of the 
curve passing through two points, but this geodesic may 
not be unique unless the two points are sufficiently close 
to one another. The problem of uniqueness now involves 
topological properties of the space V N . For example, there 
is a unique geodesic passing through two points on a 
sphere, except when the two points are at the ends of 
a diameter. In this latter ease, all great circles passing 
through the two points are geodesics. 

For a Euclidean space, referred to rectangular cartesian 
coordinates, the Christoffel symbols are zero. Hence the 
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geodesics are given by d 2 x l jds 2 = 0, whose solution is 
x i _ s B\ where A 1 and B l are constant vectors. 
That is, the geodesics are straight lines. 

§ 27. Null-geodesics 

Equations (15.3) state that 
dx i dx s 


(27.1) 


g,i ds ds 


d 1 dx i 
ds \ g,J ds ds) 


along any portion of a curve which is not null. On diffe¬ 
rentiation we obtain 

d t dxtdxA 6 / dx‘ ^ (<W\ 

ds -te-&) = ds li) ~ gii ds 6s\ds / 

It follows from (26.4) that the invariant 

is zero at all points on a geodesic. Thus the indicator e 
cannot change abruptly along a geodesic, and so if the 
tangent vector is not null at any one point, it cannot 
be null at any other point on the geodesic. On the other 
hand, if the initial direction is null, then the curve is 
null and it is of course impossible to introduce the arc- 
distance as parameter. Instead we now say that a null 
curve a?* = (£) which is a solution of the equations 

d 2 x l | il dx* dx k 

(27.2) -jp + | . fc J -fa-fo ~ 

is a null-geodesic* 

The null-geodesics in the V t with line-element (14.2) 

satisfy the equations ^ = 0. Therefore the null curves 

given by (15.2) do not satisfy these equations unless r, 
0 and y are constants. It follows that a null curve is not 
necessarily a null-geodesic. 

§ 28. Geodesic coordinates 

We shall now show that it is always possible to choose 
the coordinate system so that all the Christoffel symbols 
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are zero at a particular point. Consider a general coor¬ 
dinate system whose values at a particular point P 0 
are ^{op anc ^ introduce a new coordinate system x ( by 
the equations 

(28.1) *-«■-«„ 5)(«-—«5,). 

The index (0) attached to any entity denotes its value 
at the point f* 0 . The brackets serve to emphasise that 
this index has no tensorial significance and that the 
summation convention does not apply to it. Differentia¬ 
tion with respect to x 1 yields 


( 28 . 2 ) £-«+$„(*-*>■ 

Hence = dj. Accordingly the Jacobian deter¬ 


minant 


(0) 


is not zero, which shows us that the 


transformation (28.1) is permissible in the neighbour¬ 
hood of P 0 . On inner multiplication of (28.2) with 
dx } jdx k , we obtain 


3x*3£* 

Thus at P n 


xi dx* ii\ . dx i 

* d^ + \jn\ m {afl ~ X m ) d^- 
We differentiate this with respect to ** and obtain 
_ d*x* t i | dx n dx* , ( i \ , d 2 x* 

" r Ww a* as* + UJw (£E " 


3x*3x*' 


We now substitute in (21.4) and derive 

iv) 


That is, 
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Hence a particular system of coordinates, called geodesic 
coordinates, can always be chosen so that the Chris- 
toffel symbols arc zero at any assigned point called the 
pole. The transformation (28.1) is not the only method 
of obtaining geodesic coordinates. However, we have 
shown that one such system exists and in this particular 
system £{ 0) = 0, which means that the pole is also the 
origin of coordinates. 

The important property possessed by a geodesic coor¬ 
dinate system is this; the covariant derivatives reduce to 
the corresponding partial derivatives at the pole because 
at this point all the Christoffel symbols are zero. We men¬ 
tioned earlier in section 9 the fact that if a tensor is zero 
in one coordinate system, it is zero in every coordinate 
system, since the transformation law of tensors is linear. 
The setting up of a tensor equation often involves heavy 
algebraic manipulations. Generally, the volume of work 
is reduced by first proving the equation with respect to 
a geodesic coordinate system at its pole. It follows that 
the equation is true for all coordinate systems at this 
point. Then if this point is general, the equation is true 
at all points of the V&, 

We shall illustrate this method by proving the multi¬ 
plication law of section 24 which is satisfied by eovanant 
derivatives. Consider the tensor 


(Atf B*) rm ^ 

Let us choose a geodesic coordinate system with pole at 
P 0 , then the eovariant derivatives reduce to the familiar 
partial derivatives at P 0 . Since partial derivatives satisfy 
the multiplication law 


1 

3h 


dw dy> 


the above tensor is zero at P 0 in the geodesic coordinate 
system. Thus it is zero at P 0 in every coordinate system. 
But P 0 is a general point, therefore the above tensor is 
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zero at all points in the V N . This establishes the truth 
of the multiplication law. 


Ex. Show that at the pole P 0 of a geodesic coordinate system 


A 


i -a — 


VA t A l_ i n 
1 (£? j 


§ 29, Parallelism 

An important property of parallelism in a Euclidean 
space, which is referred to rectangular cartesian coor¬ 
dinates, is this: a parallel field of vectors A t is obtained 
throughout a Euclidean space if the components A i are 
constants. We can express this analytically either in the 
form dAJdt = 0 or by dAJBx* — 0, Since the Christoffel 
symbols are zero, we can write these equations equi¬ 
valently in the tensorial forms dAJdt — 0 or A iti — 0 
respectively. This suggests two ways of generalising the 
concept of parallelism to a Riemannian space. We shall 
see in section 31 that the partial differential equations 
A iti = 0 are in general not consistent. Thus the second 
suggested generalisation is not a profitable one. Also the 
intrinsic derivative dAJdt is only defined along a curve, 
hence, using the first suggestion, we can only define 
parallelism along a curve. Formally, the vectors A t 
constitute a field of parallel vectors along the curve 
x* = x l (t) if A t is a solution of the differential equations 


(29.1) 


dA,^dA 1 _f l \ <fe* 

6t - dt [&} 1 dt “ 


These equations form a set of A r differential equations 
of the first order, and consequently if the vector A { is 
given at any one point on the curve, it is uniquely deter¬ 
mined at all other points of the curve. We may also say 
that a field of parallel vectors is obtained from a given 
vector by parallel propagation along the curve. Since 
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dA i 

St 


“5 


U A 

di 


we can write the condition for parallelism along a curve 
in the contravariant form 


(29.2) 


dA* _dA< m dx± 
dt ~ dt + 1/fcj dt 


— 0. 


Wc see from equations (26.4) that the unit tangent 
vectors form a field of parallel vectors along a geodesic. 

The magnitude A of the vector A 1 is given by 
{ A yi- ew g ti A ( A*. On differentiation we have 


2A~ = j l (e l A)gijA i A*)=j t (e { A)gijA i A i )-2e lA) g ii AK 


This equation becomes A(dAjdt) = 0 if A‘ forms a field of 
parallel vectors. We deduce that the magnitude of all 
vectors of a field of parallel vectors is constant. 

The angle 0 between the two vectors A i and B i is given 
by ^Beosfl — g^A^BK We obtain on differentiation 

AB j t (cos 0) + B+A ^ cos B>+A‘ — ). 

When both A 1 and B‘ form fields of parallel vectors, this 
equation reduces to d(cos 0)/df = 0, provided that 
neither of the vectors A i and U 1 is null. Tliat is, the 
angle between two non-null vectors remains constant 
whilst both undergo parallel propagation along the same 

curve, * 

The vector obtained at Q by parallel propagation 
from P depends on the curve joining P to Q. Hence 
parallel propagation around a elosed curve docs not neces¬ 
sarily lead back to the initial vector. As an example con¬ 
sider the V 2 formed by the surface of a unit sphere. Un 
choosing spherical polar coordinates x l = 0, x - y> the 
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metric of the V 2 is given by ds* = dO 2 + sin 2 Odp 2 and a 
brief calculation shows that the only non-vanishing Chris- 

toffel symbols of the second kind are j } = —* sin 0 cos 0 

{ 2 ) (22J 

Til “ cot ^et us se * cc ^ ^ ie sn *all circle 0 = a for 


the purpose of parallel propagation of the vector AK 
Along this circle dOjdt = 0, and consequently the equa¬ 
tions (29*2) reduce to 


dA 1 

-cos a sin a A 2 = 0; 



cot a A 1 = 0 


whose solution can readily be obtained in the form 

A 1 — sin a [c sin cos a) + d cos cos a)]; 

A 2 — e cos (ip cos a) — d sin cos a)> 

where c and d are constants. Suppose we choose the vector 
A i to be (1, 0) at the point defined by ip = 0, Then by 
substitution we have c — 0 and d = cosee Thus the 
vector A 4 is uniquely determined by the components 
(cos [ip cos a], — sin [y> cos at]/sin a). The result of parallel 
propagation around the small circle then yields the vector 
(eos[2jtcosa],—sin [2ti cosot] /sin a) which differs from the 
original vector (I, 0), The case of the great circle a = nj 2 
is exceptional, as along it parallel propagation does lead 
back to the original vector. 

Ex* In the F* with metric ds* *=* du a 2?Mudv + dv\ where A 

is a function of u and v, show tiiat the tangent vectors to the 
curves t* = constant form a field of parallel vectors along the 
curves c — constant. 


§ 80, Covariant derivative 

By means of the concept of parallelism we are now 
able to supply a proof that the right-hand side of equation 
(23*8) constitutes a tensor. Consider a curve C determined 
by the equations = x i (t) t Choose p arbitrary contra- 


8rl£; 
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variant vector-fields X\ 1)t . • X\ T) each parallel 

along the curve C and s arbitrary covariant vector-fields 
y (1)( , y (aJ1 , .. Y {t)i which are also parallel along the curve 

C. Then 


(80.1) 

and 

(30.2) 


dt 

dY (g)i 
dt 



(fi — 1) 2, . . . p), 

(a = 1, 2,... s). 


Now let us set up the invariant 


I = ^(J) ^(2) ■ • %'&) y «*>».' ■ Y MV 

The derivative of this invariant witli respect to i is also 
ail invariant, and by application of (30.1) and (30.2) and 
appropriate changes of dummy indices we find 




dA u 


dt 


* - , [u 

, \ dz’’ 

+ £ 1 ^ ;: ^ 1 a+i 

ft) ~dt 

V ( 2 

: 1 

_ y a u I”"- , 

n] dt 


We deduce from the quotient law that the expression in 
brackets on the right-hand side of this equation is a 
tensor, which we call the intrinsic derivative, and denote 
by <M“ l "?'/&■ ^ follo ws immediately that 


( rdX££ 


£^« l .* u a -i feu a+ i 
j , r i- T i 


dt 


dt 


dxf 1 


a-1 


■-•fc) 


£ A *v ;; 'e-i tT e*i- r * (r, »)] ‘ 
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The quotient law again shows us that the expression in 
square brackets is a tensor, which we denote by A* n 
and which we call the covariant derivative* This completes 
the required proof. 


chapter v 


CURVATURE TENSOR 


§ 31 . Riemann-Christoffel tensor 

We will now investigate the commutative problem with 
respect to covariant differentiation. Let us begin with the 
covariant derivative of an arbitrary covariant vector 



A further covariant differentiation yields 


*J,np 


9m, _mw* . 

_ d&daP \jhJ dx» 1 daP \j'nj \jp\ dx n 


We interchange the indices n and p and subtract. After 
changing several dummy indices we have 

Since A $ is an arbitrary vector, it follows from the quotient 
law that the expression in square brackets is a mixed 
tensor of the fourth order, of contravariant order one and 
co variant order throe. Using the notation 



E 
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(31*1) Kjn, = ( ? J — ^ { 7 -J + { J (j-p) — {J { ? '„} * 

we observe that R l jnp is a tensor of the fourth order, 
called the Riemann-Christoffel tensor. It is formed 
exclusively from the fundamental tensor g (i and its 
derivatives up to and including the second order. This 
tensor does not depend on the choice of the vector A it 
We can now write 

(81.2) Aj t np A j tPn — Rj np Ai . 

It is clear from this equation, that the necessary and 
sufficient conditions that the covariant differentiation 
of all vectors be commutative, is that the Riemann- 
Chxistoffel tensor be identically zero. We stated in section 
29 that the equations A t j — 0 are not as a rule consistent* 
In fact, equations (31.2) show us that a necessary, but 
not sufficient condition, for the consistency of A i f ^ 0 
is that B^ fnp A l = 0, a set of equations which is not 
generally satisfied. 

Referring to the definition (31.1), we observe that 

{81.8) = -*.**»• 

That is, R l inp is skew-symmetric with respect to the 
indices n and p. 

liX. 1. 1 rove that Rj np + + R \pin ^ 

Ex. 2. Prove tliat J ? 1 lnp »* 0. 

§ 32. Curvature tensor 

We now introduce the eovariant curvature tensor 
defined by 

( 82 *D fl rtn , = gr,R.W 

On substituting from (31.1) into (32.1) we obtain 
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= [ grt i>)] _ {?>} L [ Sti U] 

+ (/I) + grt {«) U) _ gw U) (?»} ’ 

which reduces on application of (20.3) and (20.4) to 

- if £ 6 '”' rl + (D (ft.) [r ”’ ,] - 


We can reduce this further by means of (20.1) and (20.2) 
and finally obtain the important formula 


(82.2) 


, J_ / 9 2 g rP , 9 a gj» _ d 8 gm _ d 2 g jp \ 

'rjnp 2 daf 1 ‘ dx r 3x p 8x> dx v dx T dx n j 


+ s] [rp, t] — \jp, s] [rn, t ]). 

From this result we immediately deduce the relations 


(32.3) 

and 


^rlTiv = 

Rfjnv = 

l ^rlnv ~ ^nprl * 


(32.4) H r ] nv + tt rnv) + R TTi » — 0- 

We are now faced by this problem: how many distinct 
arithmetical non-vanishing components does the tensor 
R rjnv generally possess? Referring to (32,3) we sec that a 
component is zero if either r = j or n — p. Thus the com¬ 
ponents apart from sign conform to the three types 
R rlrj , R rjrv and R rlnp where r, j, n and p arc distinct from 
one another. There are as many components of the type 
R rjrj as there are ways of combining r and j. That is, 
£ A 7 (A r — i). There are as many components of the type 
R . as there are ways of combining j and p after selec¬ 
tion’of r. That is, h\ 7 (N — l)(N — 2). The number of 
combinations of r, j,n and p is -^A'(A T —1)(N— 2)(N 3). 

But R rjnv is determined, except for sign, when r is paired 
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with either n or p . Thus there arc IN(N —1 ){iV-2)(iV-3) 
components of the type R rjnv * Equations (32.3) enable 
us to rewrite (32.4) in the form 

^irnp + + Rjnpr = 

And so with a given set of four indices* only one in¬ 
dependent equation of the type (32.4) exists. Further, 
if the indices /, r, n and p are not distinct* (32.4) 
reduces to one of the equations (32.3). Hence the 
number of independent equations of the type (32.4) is 
-LN(N— l)(iV — 2 )(N — 3), Therefore the number of 
distinct components of R rJnv is 

%N(N- 1) + %N(N- 1} (N —2)+i*V{A 7 -l )(N-2)(N-S) 
- — 1)(W — 2 )(N — 8) = ^N*(N* - 1 ). 

In particular we see that the curvature tensor of a V 2 
lias only one distinct non-vanishing component. 

Ex. 1. Prove the relation (32,4) by setting up a geodesic coor¬ 
dinate system. 

d*G 

Ex. 2. Prove that /? im = — G ~— for the F, whose line- 
element is ds* = du 2 + G 2 dv 2 , where G is o function of u and v . 

§ 33. Ricci tensor - Curvature invariant 

At first sight there appear to be three different ways of 
contracting the Ricmann-Christoffel tensor R] inv * We have 
R* lnjf — S l *R S hip = because R alnv is skew-symmetric 
in $ and L We see from (31.3) that R l m = — R l nn . 
Hence we need only consider the contraction* called the 
Ricci tensor, defined by 

(33.1) R jn = R] M = g^R cint . 

On contraction of l and p in (31.1) and on substitution 
from (20.6) we find that 
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(33.2) 


R 


a» 

in dx i dx n 


o °8 

+ (.L)(;H;)s <ioeva ' 


from which it is dear that R in is symmetric. (If g is 
negative, we must replace log Vg by log V~§*) 

The curvature invariant is defined by 


(33.3) R = g /n %r 

A spaee for which R ti =* Jg^ at all points, where / is 
an invariant, is called an Einstein space. Inner multi¬ 
plication by shows that R = NJ , Thus for an Einstein 
space 

(33.4) R if = ^ Rgif 


Ex. For a V %1 prove that gR if = — &R ** — 

Hence deduce that every F a an Einstein space. 

§ 34. Blanches identity 

Let us choose a system of geodesic coordinates. On 
referring to (31,1) we have by covariant differentiation 
that 

9 i d 2 ll\ d 2 f l\ 

R inp.T = faTr (R'inp) “ (jpj dx p dx r f/rtj 

at the pole. Cyclic interchange of p and r gives us two 
other equations. We obtain by addition 

(34.1) R\i up ,.r + R l Jpr,n + ^Jrn,p ~ 

This is a tensor equation, true at the pole of a geodesic 
system of coordinates. Thus it also holds for every coor¬ 
dinate system at that pole. Further any point can be 
chosen as the pole of a geodesic coordinate system. 
Therefore equation (34.1) is true at all points of space. 
Inner multiplication by g lm yields the Bianchi identity 

(34.2) Rminp w r + ^mipr t n H" = ^ 
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The Einstein tensor is defined by 
(34,3) = — 

The inner multiplication of (34.2) by and the 

application of (33.1), (33,3), and (32,8) give us the 
equation 

R tT -g*R iT>n -g**R mr , v = 0, 
which may be written 


(84.4) R r = 2g*R jrn . 

Differentiating (34.8) eovariantly, we obtain 

That is 


(84.5) Gf fi | = 0. 

This equation is important in the theory of Relativity, 


§ 35, Riemannian curvature 

From any two vectors A { and B l at a point of a v N , 
we can construct the invariant R rinv A r A n B i B p > Let us 
consider what happens if we replace the vectors A i and 
B 1 by the two linear combinations 

X* = XA i + f*B* t Y i = pA * + t 

where X, ft, p and t are invariants. A straightforward 
calculation, with the aid of (32.3), shows that 

R rjnp X'X«YiY* - (XT- P {i)*R rinp A'A»BiBK 

Tlius the expression R rjn9 A r A n B*B*j which is an in¬ 
variant with respect to coordinate transformations, is 
almost an invariant under linear transformations of vec¬ 
tors, In order to obtain an expression which is also in¬ 
variant under linear transformations of vectors let us 
evaluate 
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(g„g»-M«)*x n r iYV D w , _ 

= {XA n + / iB n )(XA”-'rpB n )(pA v )-cB v )(pA*+TB») 

-(XAv+iiBJipA’+TBvMAfi-fiBMpAi+TBi) 

= {e A VA*+e Bf i*B* i- 2?./i cosOAB)(e A p*A 2 +e B 4 2pr cos OAB) 
-{e A XpA*+e B /nB*+[te+PPi cos 0 ABf 
= (At -fip)\e A e B - cos 2 0)A 2 B* 

= (Xz-pp)'(g Tn gi,-g r ,g in ) ArAnBIB ’’ 


where 0 is the angle between the vectors A* and B l . It 
follows that 


(85.1) 


_ R rinp A^PBrlP _ 

(firnfirt ~ grv'iiM rAnB>BV 


is an invariant which is unaltered at a point, when the 
two vectors determining it arc replaced by any linear 
combination. This invariant is called the Riemannian 
curvature of the space V"^ associated with the vectors 
A* and B'. Note that the denominator of K is unity if 
the vectors A i and B' are orthogonal unit vectors. 

At any point of a two-dimensional space there exist 
only two independent vectors. Hence the Riemannian 
curvature of a V 2 is uniquely determined at each point. 
Its value is easilv found by choosing the two vectors 
whose components are (I, 0) and (0,1} respectively. Then 


(85.2) 


£ _ -^1213 _ _ -^12 12 

fill §22 fil2 2 £ 


§ 80. Flat space 

We say that a space is flat if K = 0 at every point of 
it. From’(35.1) the necessary and sufficient condition is 

R Tinv A r A"B>Bv = 0 

for all vectors A* and B 1 . In virtue of equations (32.3) 
it follows that 


Rrjnv "T B n)rp "I" R«vri + Rrtmi — 
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That is, 

^rfnp 4 " ^Tpnj ~ 

which can be written 

^rjnp — ^tpjn * 

Interchanging /, n and p cyclically we obtain 
Thus wc have 


■®rnvi ® rinv* 


^rinp “ ^rpln ^rnpj' 


Substitution in (32,4) immediately yields 


&Ttnv — 


§ 36 


Conversely, if R rnfp = 0, then it is clear that K => Q> 
Therefore, the necessary and sufficient condition that a 
space Vn be flat is that the Ricmann-Christoffel tensor 
be identically zero. 


In a flat space R] inp — 0, hence we deduce from (31.2) 
that in a flat space the equations A ti — 0 are con¬ 
sistent, Inner multiplication by dx f jdt yields SAJdt — 0, 
Thus in a flat space the property of parallelism is in- 
dependent of the choice of a curve* We may therefore 
say that parallelism is an absolute property of a flat 
space* 

A familiar example of a flat space is the Euclidean 
plane for which the metric is = dx 2 + dy 2 in rec¬ 
tangular cartesian and ds 2 = dr 2 + r 2 d0 2 in polar coor¬ 
dinates* 


Ex* If the metric of a two-di mens tonal flat space is /{^[(dx 1 ) 1 + 
(<&’*)*]> where (r)* = (**)> + (x*)\ show that /(r) = c(r)\ 
where c and k arc constants* 


§ 37, Space of constant curvature 
Let us now investigate spaces in which the Riemannian 
curvature at every point does not depend on the choice 
of the associated vectors A* and B ** From (35.1) the 
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necessary and sufficient condition is that 

{K( SrnSip - g rp g in ) - R rtnv }ArAWB* - 0 

for all vectors and BK A similar calculation to that 
of the previous section will show that this condition 
reduces to 

^r}np 

where K is now a function of the coordinates x l t 
Covariant differentiation gives us 

RrinV't “ ^ t t(Sm8ip Srplbn)* 

We substitute this result in Bianchi’s identity (34*2) and 
obtain 

^ p r(^mn Bjp BmpBin) 4"* ^ ,n(BmpB}r BmrBjp) 

4“ -^ r j>(gmrSin BmnSjr) 

Inner multiplication by g mn g fp yields 

(N — l)(N — 2 )K ft = 0* 

Hence if N > 2, it follows that K is constant. Thus we 
have proved Schuris theorem; 'if at each point of a space 
Vn, (N > 2), the Riemannian curvature is a function of 
the coordinates only, then it is constant throughout the 
IV* Such a V N is called a space of constant curvature. 

The metric of the V % formed by the surface of a sphere 
of radius a is ds 2 = a 2 (d0 2 + sin 2 8dy) 2 ) in spherical 
polar coordinates. The reader is asked to verify that 
I£ ma = a 2 sin 2 0 and that it follows from (35*2) that the 
surface of a sphere is a surface of constant curvature 1/a 2 . 

Ex, 1. Show that a space of constant curvature is an Einstein 
space. 

Ex, 2, In a Euclidean V tt prove that the hypersphere 
x l — c sin 0 sin tp sin rp % ar* = c sin 0 sin tp cos y 
x 3 — c sin Q cos sc* = c cos 6 

is a of constant curvature 1/e*. 












CHAPTER VI 


EUCLIDEAN THREE-DIMENSIONAL 
DIFFERENTIAL GEOMETRY 

Euclidean Geometry investigates the properties of 
figures which are invariant with respect to translations 
and rotations in space* It may be subdivided into Al¬ 
gebraic and Differential Geometry, The former studies by 
algebraic methods the theory applicable to entire con¬ 
figurations such as the class or degree of a curve. The 
latter discusses by means of the Calculus those properties 
which depend on a restricted portion of the figure. For 
example, the total curvature of a surface at a point only 
depends on the shape of the surface at that point. Suc¬ 
cinctly we may say that Differential Geometry is the 
study of geometry in the small . This chapter is not intend¬ 
ed to be a complete course on the subject. However, 
sufficient theory is developed to indicate the scope and 
power of the tensor method* 

§ 88. Permutation tensors 

In the Euclidean three-dimensional space let us in¬ 
troduce the quantities defined by 

( 38 * 1 ) £ iik = y/g e m \ e m = —— 

where e iSk arc the permutation symbols defined by the 
following conditions;- 

(i) e iik — 0 if any two of the indices i, j and& are equal, 

(ii) e U3 — e 2 $i ~ e ziz “ + L 

(iii) e 13a = #32i = — — 1 , 
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and g is the determinant formed from the fundamental 
tensor g fj of the space referred to some general coordinate 
system, which is not necessarily rectangular cartesian. 
The definitions show us that e m , e m and are skew- 
symmetric in all their indices, (Note that in this chapter 
the range of Latin indices is now from I to 3), 

We shall first prove that although e ijh is not a tensor, 
both s m and s iik are tensors. We observe that 

dx l dx* dx* dx* dx k dx* dx* dx k 

c, ' , * ai 1 ai™ ~ e ‘ iK Ci,k d£™W i d^' 


dx^ dx k # * 

Thus e ijk —7 is skew-symmetric in l and m , Sinii- 

larly it is skew-symmetric in the indices I, m and n* But 
this expression apart from sign is the Jacobian determi- 
dx r 

. It therefore follows from the theory of deter- 


mint ■=— 
ox* 

minants that 


3 #* dx* dx k _ ]3# r 

e<ili dx™ dF^ ~ Clmn \W 


Now the fundamental tensor g iS transforms to g tj wiien 
we transform to the coordinate system xK We find from 
( 21 * 1 ) that their determinants satisfy the equation 


I = g 


3af * 
dx* 








Hence 


Re 1, 


dx* dx* dx k 
= fa™ 3i" 


3a? 
dx* dx j dx k 
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from which we see that e m is a eovariant tensor of the 
third order* Also 


#imn _ 

" VI 

therefore, 




1 \d& 


i Sx 1 dx m dx n # 

e ‘ mn = VI e<ik d&d&d ** 1 


eiik 


dx l dx m dx n 
dx l dx* dx k 1 


which shows that 8 tmn is a contravarianfc tensor of the 
third order. We call a iik and e m the permutation 
tensors. 

If the coordinate system is rectangular cartesian, g — I 
and the permutation tensors have as their components 
the permutation symbols* In this coordinate system the 
covariant derivatives s iik z and a iik which are both 
tensors, are zero. Thus the covariant derivatives s ifk f 
and a iik ( are zero tensors in all coordinate systems. 
Hence the permutation tensors behave like constants 
with respect to covariant differentiation. 

From a covariant vector A g we can form the vector 

& = 

and we call B* the curl of the vector A i and write it 
curl A P 


Ex, 1* Prove directly that — 0. 

Ex. 2, Prove that £ iik - 

Ex. 3. Show that the components of curl A ( are 



W.\ 

1 

/i/Jl 


1 

(M* 



J ’ 

VI 

\lx s 

7>x' ) ’ 

Vg 

\lx' 

_ 

i.r s } 


§ 89, Vector product 

We can form the contra variant vector 
(39.1) C* = 4*AsB k 

from the two covariant vectors A { and Bp To find the 
geometrical interpretation of the vector C*, let us choose 
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a rectangular cartesian coordinate system. In such a 
system the distinction between contravariance and 
covariance disappears and the C* are now the components 
of the vector product* of A t and Bp Thus C* is a vector of 
magnitude ^Usinfl, orthogonal to both vectors A £ and 
B it where 6 is the angle between these vectors. Its direc¬ 
tion is uniquely determined by the fact that A it B t and 
C H form a right-handed system. 


§ 40. Frenet formulae 

In this section we shall investigate the theory of 
twisted curves. Let the curve be given by the equations 
x l — #*($), where the parameter s measures the arc- 
distance along the curve. Then the unit tangent vector 
T 1 to the curve is T* = (Itf/ds and satisfies the equation 

ST* 

g t/F'T* = 1 . On differentiation we obtain g^-^-T* — °* 

That is, dT l l&$ is a vector orthogonal to the tangent 
vector. Its magnitude k is called the curvature of the 
curve, and the unit vector 


(40.1) 


»‘= 1 T 

K 08 


is called the unit principal normal vector. We define 
the unit binormal vector B i to be the unit vector or¬ 
thogonal to both the tangent and principal normal vec¬ 
tors, and oriented so that the tangent, principal normal 
and binormal vectors form a right-handed system. Hence 
from (39,1) we have 

(40.2) B i = B iik T j N k . 

Since the vectors T l and N* arc orthogonal, g^TW* = 0 , 
We differentiate intrinsically and obtain 


<5A 7 * 6T l 


* D* E. Rutherford, Vector Methods, p, 7, 
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dT* 

We substitute for from (40.1) and then replace 
ds 

gby both being equal to unity, and obtain 

(8N* 




Now we differentiate = 1, which gives us 

gijN* -y- = 0, Thus 


«»"■ (f +- °- 


gift 

Therefore the vector -r-1- kT i is orthogonal to both 

the tangent vector T* and the principal normal vector 
NK Accordingly it is in the direction of the binormal 
vector IP, and we can write 

<‘°' s > B ‘ = 7 (iT + * TI )' 

The invariant x thus introduced * is called the torsion of 
the curve. Note that it may be positive or negative in 
contrast to the curvature which is essentially positive, 
being the magnitude of a vector. 

We differentiate (40,2) intrinsically and substitute 
from (40.1) and (40.3) and the result is 

52_ e »^iw, + c </.r 

ds ds o$ 

= KeWNjNt + E^T $ [tB k — *T k ] 
which reduces on account of the skew-symmetry of e ijk 


to 




But T if N t and Zf ( form a right-handed system of unit 
* Some text books use x to denote the reciprocal of the torsion. 
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vectors; and so s iik Tj B k *= — N* and our equations become 
(40.4) tV<. 


Equations (40.1), (40.3) and (40.4) are called the 
Frenet formulae. On account of their importance in the 
theory of curves wc group these formulae together for 
convenient reference in the form 


(40.5) 


dr* 

ds 


dN* 

ds 


= — kT' 


KN* 


+ xB‘ 


SB* 

~ds 


— T N*. 


If the coordinate system is rectangular cartesian, the 
intrinsic derivatives become the ordinary derivatives and 
we have the well-known Frenet formulae. 


Ex. I. Show that «* = g it ——. 

5s (5s 

Ex. 2. Prove that t = N, — . 

ds 

Ex. 8. A helix is defined to be a curve whose tangent vector makes 
a constant angle with a fixed direction. Prove that the 
necessary and sufficient condition that a curve be a helix 
is that the ratio of the torsion to the curvature is constant. 

Ex, 4. Prove that & 1ft — 1 ? = * 6 — ^. Hence a curve 

ds ds 4 ds 3 da \te J 

is a helix if and only if — ---- — 0. 

ds ds 1 6 s* 


§ 41, Surface - First fundamental form 
We will now investigate the theory of surfaces. The 
three equations 

(41.1) w l — x l (u l f u 2 ); x 2 = x*(u\ u 2 ); a? — x*(u l , u 2 ). 
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where u l and u 2 are parameters and the x i are three func¬ 
tions of u 1 and u 2 which are real and continuous, generally 
represent a surface. These equations may be briefly 
written m* = x*(u a ) if we adopt the convention that Greek 
indices will always have the range 1, 2, A point at which 


the Jacobian matrix 


,x M 1 


is of rank two is called a regular 


point* A point may fail to be regular either because it 
is a singularity on the surface, for example, the vertex 
of a cone, or because it is a singularity of the parametric 
representation, for example, the poles of a sphere, (See 
equation (41.2))* 

Each pair of values of u a determines a point on the 
surface. That is, the u a form a coordinate system upon it* 
Thus an equation of the type /(t* 1 , u 2 ) = 0 must define a 
curve. Although a unique point on the surface corresponds 
to a fixed pair of values of the converse is not neces¬ 
sarily true. This is illustrated by the equations 


(41.2) x l — a sin^ 1 cosu 2 ; x 2, — a sin u l sin u 2 ; x 2 = a eosu 1 , 


which specify the surface of a sphere of radius a. In 
practice we restrict the range of the parameters so that 
the converse will also be valid. In our example, the 
parameters are limited to the ranges 0 <Ztz and 
0 £u 2 <2tt. Then, except for the two poles, there 
corresponds to each point of the sphere a unique pair of 
values of u a . The equations u l = 0 and u l — n represent, 
not curves, but the poles respectively. At the poles, the 
coordinate u % is indeterminate. 

We shall exclude singular points from our discussion 
by considering only that portion of a surface on which 
there is a unique correspondence between its points and 
pairs of values of the coordinates u a . That is, in this 
portion every curve of the family u l = constant intersects 
every curve of the family u 2 — constant in one point 
only. The curves u l = constant arc called u 2 -curves 
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and the curves u % = constant the a 1 -curves. Collec¬ 
tively they are designated the coordinate or parametric 
curves on the surface. Along a coordinate curve we 
choose the positive direction to correspond with increasing 
values of the variables u 1 or u 2 respectively. 

The contravariant vectors dx* and du a which represent 
the same displacement in space and on the surface respec¬ 
tively are connected by the equations 

9 * 27 ^ 

(41,8) dx* = du a 

where we extend the summation convention to apply to 
Greek indices. Hence the line-element ds on the surface 
is given by 

ds 2 = g ti dx'dx> = g« ~ ~ du* du/>. 

Let us write 

, . dx* dx 1 

(41.4) = 


from which it is clear that a a $ is symmetric, and so we 
have 


(41.5) ds 2 = a a fidu a dufl' 

To this equation we apply the quotient law. Since a ai $ 
is symmetric it follows that n a fi is a eovariant tensor with 
respect to transformations of the coordinate system u®. 
We call it the fundamental surface tensor. Also 
a^du^ditfi is named the first fundamental form of the 
surface. 

Let us select rectangular cartesian coordinates in space. 
Then we readily obtain the metric in the well-known form 

ds 2 = Edu 2 + 2 Fdudv + Gdo\ 


where u = u\ v — u 2 and 


dx* dx* _ ^ dx* dx* p ^ y ^ 

du du * ^ du dv * i dv dv 


F 
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§ 42. Surface vectors 

When we transform coordinates in space, das* is a 
contravariant vector but du* is an invariant. On the other 
hand, if we transform the coordinates u a on the surface, 
dx i is an invariant but du A is a contravariant vector. 
Thus the equations (41.8) indicate that we can regard 
dx i jdu a as both a contravariant space vector and a co¬ 
variant surface vector. We may therefore introduce the 
notation 


(42.1) 


04J* 


and rewrite (41.4) in the form 


(42.2) a af = 


A curve on the surface is represented parametrically 
by the equations u a = w®(£). The vector du a jdi is a tangent 
vector to this curve. Its space components arc then given 
by the equations 


(42.3) 


dx i Star* du a 
dt du a di 


du a 


But if the components of the space tangent vector dx i jdt 
are fixed, (42.3} consists of three equations for the two 
unknowns du a /dL They are not consistent unless the vee- 
tor lies in the surface, when a unique solution would exist. 

Next we consider a surface vector-field A a . We can 
set up a unique curve C on the surface by the differential 
equations du a /dt = provided that the vector-field is 
fixed at some particular point. Then A a is the surface 
tangent vector to this curve C, Let us designate the space 
components of A a by A\ and so equations (42.3) state 
that these components are connected by the relations 

A 1 — A** 


(42.4) 
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The magnitude of the vector A { is given by 

M»* - g 

That is, 

(42.5) (A)* = a ap A*A$. 

In particular, du a jds is the unit tangent vector to a curve 
on the surface if the parameter s measures arc-distance 
along it. 

The angle 0 between the two unit vectors A { and B l 
is obtained from 

cos Q = g ii A i B i = g (f x { a x^A a Bf. 

That is, 

(42.6) cos 0 = a a fiA a B0> 

It follows that the necessary and sufficient condition for 
the orthogonality of two surface vectors A a and B& is 

(42.7) a A fiA*Bfi = Q. 

We sec from equations (42.5), (42.6) and (42,7) that 
the familiar formulae apply equally well on the surface 
provided that we employ in them the components of the 
vector in the surface and the surface fundamental tensor. 
Also we can raise and lower indices of surface tensors in 
the usual way with the surface fundamental tensor a a $ 
and its conjugate symmetric tensor a a $. These two tensors 
are connected by the equations 

(42.8) a a/l a*v = gr, 

where d is the two-dimensional Kronecker delta. It is 
worthy of note that 

(42.9) a 11 — a 2i la; a 13 — a 21 — — a 12 /a; a 22 — 
where 

a — a 22 (%)*• 

Ex, Prove that A a = z a r A t , 
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§ 43. Permutation surface tensor 

In section 38 we introduced the permutation tensors 
in space* Similarly we introduce on the surface the 
quantities defined by 

(43.1) e aff = */a S <' a fi, 

where 

^11 ~ ^22 = ^12 ” H “ ^21 = 

It is left as an exercise for the reader to show that e&p 
and are a co variant and a contra variant skew-symme¬ 

tric surface tensor respectively. They are called the 
surface permutation tensors. We see that s a p can 
be obtained from by lowering the indices since 
— a ar a 0t£y*. 

Now wc derive an important formula for the angle 0 
between two unit vectors A a and B*. This angle is given 
by cos 0 = a a pA a BP* Accordingly, using (42.6) 

sin 2 0 = 1 — a a flA a BPctyjAy B* 

= — aa/tayt)A*AYB*B S 

= ae^ e y fl A* A? B*B S 

= eeieyfiAvArBfiB* 

= 

To remove the ambiguity in sign the convention is made 
that we choose that value of 0 which satisfies 

(43.2) sin 0 = + e ad A a B s . 

In accordance with the convention just made, we say 
that the rotation from O to i> is positive if the invariant 
E a fiC a D# is positive. This in effect chooses the positive 
rotation as that one which rotates O to D a through an 
angle less than or equal to n . 
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Let us form the contravariant vector 
(43.8) B# = && A a 

from the covariant unit vector A a . Its magnitude is 
given by 

(U)* = a a fiB*Bt - a m0 ^f^A r A i 

= + 2a^e'H*'A 1 A 2 + a ll {e*mA 2 ) 2 

— “ 0*22 Ml) 2 + ^llM*)*}* 

That is, in virtue of (42.0) 

M) a - cP(A t )* + 2a 12 A l A 2 + a**(A t )* 

= a a #A a Ap « (^ 4) 2 — 1 , 

Thus B a is a unit vector. Further the angle 0 between 
A a and B* satisfies 


sin 0 = SapA* Bfi = e a ^A^ Bp = B& B$ = I, 


Therefore 0 = sr /2 and by the above convention, equation 
( 43 . 8 ) determines the unit vector B a orthogonal to the 
unit vector A a and oriented so that the rotation from 
A a to B a is positive. 

We now apply (43.2) to calculate the angle m between 
the coordinate curves. The unit tangent vectors to the 

u 1 - and u a -curves are ~= Sf and *_ respectively, 

V<hi v%2 

Hence the angle m satisfies 


(43.4) sinct> — 


V<hi& 




22 


y/a n a. 


H2 


22 


-V—. 


We also note from this equation that the rotation from 
the direction of a u l - curve to a u 2 -curve is always positive. 
It is easy to deduce that the necessary and sufficient 
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condition that the coordinate curves are orthogonal 
everywhere on the surface is that a 12 vanishes everywhere. 
Wc then say that the coordinates are orthogonal 
curvilinear* 


g 44. Surface covariant differentiation 
We can form Christoffel symbols starting with the 
fundamental surface tensor and hence introduce 
covariant and intrinsic derivatives which are now sur¬ 
face tensors. So $ will denote the covariant derivative 
of A a with respect to vfi. Written in full 



There will be no confusion between the Christoffel sym¬ 
bols formed with the g if and those formed with the tensor 
a*# since the Latin and Greek indices will clearly dis¬ 
tinguish which symbol is meant. We can also form the 
Riemann-Christoffel surface tensor Ji^ ^ from and the 
curvature surface tensor R € p r $ = R a $ y $ , 

Following the argument of section 26, we can show 
that the geodesics on the surface are the solutions of the 
differential equations 


(44*1) 


d fdu*\ _ d 2 u a ( | AvP duv 

fo \dT) = ~d^ " r [ffy) ~ds ~ds 


Corresponding to (27.1) the geodesics also satisfy the 
equation 


(44.2) 


du a du$ 

a ' afi Is ~ds ^ 1 ‘ 


Thus, in practice, we need only consider one of the two 
equations (44*1) together with the first order equation 
(44.2). 

As in section 28 we can introduce a system of geodesic 
coordinates on the surface so that at any particular point. 
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called the pole, all the surface Christoffel symbols are 
zero. At the pole, covariant and intrinsic derivatives 
reduce respectively to the corresponding partial and total 
derivatives* 

Again the theory of parallelism, outlined in section 29, 
applies to surface vectors. The vector-field A* is said 
to be parallel along the curve u a = u*(t) if 


dA* _ dA* 
~&t - dt + 



= 0 * 


Note that the covariant derivatives of a a $ and 
are zero. We cannot apply the method of section 68 to 
prove that s a p t7 and £*$ are both zero because it is general¬ 
ly impossible to choose a rectangular cartesian system of 
coordinates on an arbitrary surface. Instead we select a 
geodesic coordinate system* At its pole, the Christoffel 
symbols are zero. Hence the partial derivatives of a a p 
are also zero. Thus the partial derivatives of the deter¬ 
minant a are zero at the pole. That is and are 
both zero there. It follows immediately that these tensors 
are zero at every point on the surface and in all coor¬ 
dinate systems. Hence the permutation surface tensors 
behave like constants with respect to surface covariant 
and intrinsic differentiation. 


Ex. Show that the conditions that the tt 1 -curves and the u^curves 
be geodesies arc j^J = 0 and j^J = 0 respectively. Simplify 

these conditions if the coordinate system is orthogonal cur¬ 
vilinear. 


§ 45. Geodesic curvature 

We discussed in section 40 the theory of the twisted 
curve in space. To this we now add the corresponding 
theory applicable to twisted curves lying on surfaces. 
Let the curve be given by the equations = u*(#), 
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where s measures are-distance along the curve. The unit 
surface tangent vector is 


(45.1) 



Since i a is a unit vector, we have a a ^t a 0 — 1 and intrinsic 

differentiation yields us = O.This shows that dtPjds 

is a vector on the surface orthogonal to the tangent vector 
f*. Let us denote the unit vector in the direction of <5 t*!6$ 
by n*. Then 

Si* 

(45*2) ~ = <m* 

OS 


where a is an invariant called the geodesic curvature 
of the curve. We also call n a the unit normal vector 
in the surface to the curve, and choose its direction so 
that the rotation from t a to n a is positive* Thus 


£ a p t a n* = 1, 


and equations (45.2) determine a uniquely in sign as 
well as magnitude. It is now clear from (43.8) and sec¬ 
tion 17 that 


and 


nP = + ^ ta 


t$ — — £ a #n a . 


Intrinsic differentiation of the first equation gives us 
dnfi mS 6t a - 

__ = — — dir w 

6$ ds 


Combining this with (45,2) we have the surface Frenet 
equations of a curve 


(45,3) 


<5f* 

6s 


= an* 


6n a 

17 




Although t a is the same vector as T l , their components 
being connected by the relation T l it is important 

to note that n a is in general neither the principal normal 
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N* nor the binormal B*. It does lie however in the normal 
plane of the curve determined by N* and B*. 

Along a geodesic of the surface dtPjds — 0 and hence 
& — 0 * Conversely if a = 0 , then 6t a j6s — 0 and so the 
curve is a geodesic* Therefore, the necessary and sufficient 
condition that a curve be a geodesic is that the geodesic 
curvature be zero. 

Ex. Prove that the geodesic curvatures of the coordinate curves are 

§ 46. Normal vector 

We shall now derive an expression for the unit normal 
vector at any point on a surface. Let us choose its 
orientation so that the i^-curve, the u s -curve and the 
normal at the point form a right-handed system. The 
surface unit tangent vectors to the coordinate curves are 
1 l 

— — 6f and ■ - 6t respectively, and the corresponding 

V a n V^aa j ^ 

space components are _ and - *— #£$2 * 

1 i . v a u 

—=.#* and In virtue of (43.4) and the co- 

2 

variant form of (89*1) we have 

1 / —-— — Efi* - , 

f &H&22 V a ll a 22 

That is, 

(40.1) I* = 

It is not clear that is a covariant vector because y/a 
appears in this equation. However the vector form is 
clearly seen from the equivalent equation 

(46.2) f, = 
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For purposes of calculation (46,1) is the more suitable 
and can be readily written as a determinant, but (46,2) 
is preferable for theoretical purposes. Since the vector 
does not lie in the surface, there is no corresponding surface 
vector We shall often require the important equation 

(46.8) gaP 4 = 0, 

which expresses the fact that the normal f* is orthogonal 
to the surface vector a^. This equation is also an im¬ 
mediate deduction from (46*1), 

Ex* The tangent surface of a curve is defined to be the surface 
generated by all its tangent lines, called the rulings* Show that 
the normals to the surface along a ruling are all parallel to the 
binomial of the curve at the point of contact of the ruling* 

| 47* Tensor derivatives of tensors 
In the theory of surfaces we require tensors which 
possess both Latin and Greek indices, for example x 4 a . All 
such tensors in this chapter will be contravariant with 
respect to space but covariant with respect to the surface* 
We can then select A* a as a typical tensor* When we 
change the coordinate systems both in space and on the 
surface the transformation law is 

it - A i?*du> 

“dx* du a ' 

We now ask what tensors can be constructed by differen¬ 
tiation? Here we follow section 30 closely* As our space is 
Euclidean the concept of parallelism in it does not depend 
on the selection of a curve* Accordingly choose an ar¬ 
bitrary parallel vector-field X i in space. On the surface 
take an arbitrary parallel vector-field F« along the curve 
C whose parameter is t Then 
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and 



Form the invariant A* a X t Y a , By differentiation and the 
use of the equations of parallelism we have 



+ 



* 

__ Va 

di 


— A*i 




That is T on making an appropriate change of dummy 
indices, 


d 

dt 


[AiX.Y*] 





Applying the quotient law wc see that the expression in 
square brackets is a tensor, which we call the tensor 
derivative of /!£ with respect to uP and we denote it 
by the semi-colon notation 

(«-i) 

It is possible to choose a rectangular cartesian system 
in space and a geodesic system on the surface. Then at 
the pole, the tensor derivatives are the partial derivatives* 
Consequently the laws of tensor differentiation are the 
same as those that apply to covariant derivatives* 

We must now extend the concept of tensor differ¬ 
entiation further to space tensors and to surface 
tensors* It is clear that the tensor derivatives of surface 
tensors are identical with their covariant derivatives* 
Following the method of this section, we see that the 
tensor derivative of a space tensor with respect to u a is 
the tensor obtained by inner multiplication of its 
covariant derivative with respect to x l by the tensor x l a . 
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For example, = A if k x&. Thus the tensor deriva¬ 
tives of g i} > g**, <5J, e iSky 6 ^* a a 0 > a*P, <5J, £ afi and are 
all zero. That is, they may be treated as constants with 
respect to tensor differentiation* 


§ 48. Second fundamental form 


By 

(48*1) 


tensor differentiation 
j d 2 x* 

Xa '*0 du a duP 



which, shows that is symmetric in a and /}. That is, 

~ x fca' Now tensor derivation of (42*2) yields 


&A ,A + gtAAv ■ °- 

We subtract this equation from the sum of the two 
similar equations obtained by cyclic interchange of a, (i 
and y. Because x\,,j is symmetric the result is 

( 48 - 2 ) = °* 

This shows that is a contra variant space vector 
orthogonal to all vectors x? r lying on the surface. Thus 
it is co-directional with the normal vector f*. Hence 
quantities b a ji must exist so that 

(4S-8) = 

Further, it follows that b a # form the components of a 
covariant symmetric surface tensor. Equations (48.8) 
are known as Gauss's formulae. Since is a unit 
vector, inner multiplication of (48.3) by and sub¬ 
stitution from (46.2) yields 


(48.4) b a $ * e ijk ^x? y .r* — Biik ^^" 

The quadratic form 

(48.5) 


du* dufi 
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is called the second fundamental form of the sur¬ 
face, It is now possible to construct the invariant 

(48,6) H = ba*0b a/ !, 

which is called the mean curvature of the surface. 

Ex. If the space coordinates are rectangular Cartesian, show 
that 

i a*»‘ , * 

55*551 

g 49, Third fundamental form 

Tensor differentiation of the identity g* = 1 yields 
us = <h That is, is a contravariant space 

vector orthogonal to the normal vector. Accordingly it 
lies in the surface. Hence quantities tj& exist so that 

(49.1) £% - 

The quotient law then states that 7$ form a mixed surface 
tensor. We now differentiate (46.3) tensorially and obtain 

which reduces on substitution from (49.1) and (48.3) to 

gtiVlrfyrft + g = 0 * 

We apply (42.2) and the result is 

(49.2) b aff = —a fiy 7jl. 

Inner multiplication by aP* yields 

(49.3) rf a = —a?*b^i 

and so we can rewrite (49.1) in the form 

(49.4) — 

These equations arc known as Weingarten's formulae. 
Introduce the symmetric surface tensor 

(49.5) C afi = g„£f 
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and form from it the quadratic form c a pdw*du$ which is 
then called the third fundamental form of the surface* 

Ex. Prove that = GY*b a y . 


§ 50* Gauss-Codazzi equations 
Wc are now in a position to obtain the central formulae 
in the theory of surfaces. For the moment, choose the 
space coordinates to be rectangular cartesian and the 
surface coordinates to be geodesic* Then at the pole, 
tensor differentiation of (48,1) yields 


,( 9 3 a* 

9 f a \ { 

a '0 r dupdutdu? 

duy U" 1 ' 

Thus 


* < r ^ h 

M_ L Ml 


VI dur (a/ijj 


The expression in square brackets is the surface Ricmann- 
Christoffel tensor at the pole. Therefore we have set 
up a tensor equation 


(50.1) 




X *\y$ — 


which must be true at every point of the surface and in 
all coordinate systems. 

Wc substitute from (49.4) in the tensor derivative of 

(48.3) and obtain 


X z;0y ~ y £* b a 0 b yc X ^* 

We can consequently write (50,1) in the form 


— b «yifi)¥ — ^(baflby* ~ KyhMr = ^aflyK * 

Inner multiplication by f* and gijX* p in virtue of (42.2) 
and (46.3) yield tiie respective equations 

(50.2) b afi:r — b a7;fi = 0 

and 

(50*3) Hpxfly = bgybpp — b a fibpy* 
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Tiie reader is asked to verify that (50*2) consists of 
only two independent partial differential equations* They 
are called the Codazzi equations. As there is only one 
distinct component of the curvature tensor in two dimen¬ 
sions, (50*3) reduces to the single equation 

(50.4) &i 2 i 2 — b±ib 2 2 ) 2 

which is called the Gauss equation. By means of (35.2) 
we can write this equation 

(50.5) K — 

where b is the determinant formed from the b a $ and K 
is the Rtemannian curvature of the surface. On a surface 
it is more usual to call K the Gaussian or total cur¬ 
vature. 

It can be proved* that a surface is uniquely deter¬ 
mined except for a translation or rotation in space when 
the first and second fundamental forms are given* This 
theorem can be precisely formulated as follows;- If a a # 
and b &0 are given functions of u 1 and u 2 , then there exists a 
surface x 4 = s?*(u a ), uniquely determined except for its 
position in space, which has a a $ du a du$ and b a $du a duP as 
its first and second fundamental forms respectively, 
provided that a^dit^d^ is positive-definite and that 
and b a $ satisfy the Gauss-Codazzi equations* 

Ex. Prove that c a 0 = 2IIb a $ — Ka a fi . 


§ 51. Normal curvature-asymptotic lines 

On a surface consider the curve u a = u ff (s), where s 
measures arc-distance. The equations of the curve in 
space will be x l — Then the space vectors TN* 
and B* and the surface vectors t* and n a at any point of 
the curve satisfy the Frenct formulae (40.5) and (45.3). 

* L. P. Eisenhart, Differentia! Geometry, p. 157—159* 
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The tangent vectors T' and i a arc connected by T* = t?xi. 
Intrinsic differentiation yields in virtue of (48.3) 


^H_ = T i ^. = T i x 1 — = T‘ — 

As 0 fa 10 d* 


ds 


> du0 . . t duP 

= $<-&+*<*-& 


Applying the Frenet formulae we have 
kN* = aw* 4 + b ai) i*tfiP, 

which becomes 


(51.1) kN* = an 1 + t* IP? 


when we designate the space components of n a by n f . 
Let us introduce the angle 0 between the principal normal 
N* and the surface normal f f * Then inner multiplication 
of (51.1) by ii yields 

(51*2) k cos 0 = b aff t*t0 t 

since is orthogonal to the vector n* which lies in the 
surface. The invariant b a $t a tP is the same for all curves 
which have the same tangent vector i° at the point on 
the surface* Accordingly wc deduce Meusuier’s theorem 
'For all curves on a surface which have the same tangent 
vector, the quantity k cos 0 is cons tan t\ This quantity 
is called the normal curvature at the point and is 
denoted by K( n }* Hence 


(51.8) 


K {n) = b afi t* t0 — 


b a fi du a duP 
a a fi du a duP' 


If we choose a plane section through the normal to the 
surface, then 0 = 0 or jf. That is K( rt ) — ± k * So the 
normal curvature in any direction is equal in magnitude 
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to the curvature of the normal plane section of the surface 
in that direction* 

Along a geodesic a — 0, therefore equations (51,1) 
reduce to kN { — bafifctPS*. Thus either k — 0 oxN* — 

We deduce that a geodesic on a surface is either a straight 
line or is a curve whose principal norma! is co-directional 
with the surface normal at every point. Conversely, if 
A M — ± ft then inner multiplication of (51*1) by n f 
yields a = 0* That is, the curve is a geodesic. 

The directions at a point on a surface which satisfy the 
equation b a fidu a duP — 0 are called the asymptotic 
directions* If all the points of a curve have asymptotic 
tangent directions, the curve is called an asymptotic 
line* The asymptotic lines on a surface are given by 
b a fidu a du0 = 0. Thus along an asymptotic line wc have 
kN £ = an*. Hence, since JV f and n* are both unit vectors, 
either k = a = 0 or the curvature and geodesic curvature 
of an asymptotic line are equal in magnitude and its 
principal normal lies in the surface. Consequently the 
binomial of an asymptotic line, which is not a straight 
line, is co-dircctional with the surface normal. The con¬ 
verse is also valid. 

Ex* Prove Enneper's formula that the torsion r of an asymptotic 
line is ^ v — K where K is the Gaussian curvature of the 
surface* 

§ 52, Principal curvatures - lines of curvature 

The normal curvature #q„) of a surface in the direction 
t a is given by *: {n ) = b a fit*tP, where the tangent vector 
satisfies a a fit*tP = 1* The maximum and minimum values 
of K( n ) can then be determined by the method used in 
section 19. They correspond to the principal directions 
determined by b a fi and arc given by the roots of the deter- 
minantal equation | b a fi — Xa afi \ — 0, which reduces in 
virtue of (50*5) and (48*6) to 

(52.1) A 2 —2 m + K = Q. 


G 
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The roots k ( ij and of this equation are called the 
principal curvatures of the surface at the point. The 
principal directions and corresponding to the prin¬ 
cipal curvatures at the point respectively satisfy 

(baft a ap)t% 1 ) — 

[b a ff ^(2) 

A point at which is called an umbilic* At 

all other points section 19 tells us that and l* 2) are 
orthogonal to one another. A curve on the surface which 
is a principal direction at all of its points is called a 
line of curvature. 

At an umbilic, equation (52.1) has coincident roots. 
That is II 2 = K and the reader may verify* that this 
result can be written 

4a(& n b 12 2 )“ 

“T [^ix( fl n^22 ^u) ^laf 0 !! ^i2~ a ia^ii)] 2 = 

Since a a fidu a du$ is positive-definite, a is positive, and we 
deduce tiiat 

°uAl2 = ^11^22 ^22^11 “ 

and so 

bn _ ^12 _ ^22 * 
a n a l 2 a 22 

Equation (51,8) now shows that K {n ) is independent of 
the direction du a jds . That is, at an umbilic, the normal 
curvature is the same in every direction**. 

Ex, 1, Prove that the lines of curvature on a surface are given by 
$y* a ay du a dtt# = 0 , 

Ex. 2, If the coordinate curves are lines of curvature, prove that 
flit = 6 h = 0, and conversely. 

Ex. 3, Prove that a surface, ah of whose points are umbilics, is a 
sphere or a plane. 

* I,. P. Eisenhart, Differential Geometry, p. 119. 

** The reader may find it interesting to relate some of the 
results of tins chapter to chapter II of Rutherford’s Vector Methods. 
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CARTESIAN TENSORS - ELASTICITY 

§ 53. Orthogonal transformations 
In this chapter our aim is to present the Theory of 
Elasticity and so we restrict ourselves to a Euclidean 
space of three dimensions. In it we choose a right-handed 
system of rectangular cartesian coordinates and we denote 
them by y it where Latin indices have the range 1 to 8. 
The line-element ds is then given by 

(53.1) ds 2 = dy i dy { = & u dy i dij f , 

where da is the Kronecker delta. 

The linear equations, (compare (6,3)), 

(53.2) y t = a ij y s + b it 

where b t form a set of three constants and a set of 
nine constants define a transformation to a new coor¬ 
dinate system y t * The necessary and sufficient conditions 
that the y t form a set of rectangular cartesian coor¬ 
dinates is 

ds 2 - dy,dyt = a^a^dy^y^ = & jlc dy,dy k . 

That is, 

i a a a a —" &sk)dtf} dlfk — 0 
for all values of dy f . Hence 

(53.3) a u a ik — $ik* 

Inner multiplication of (53,2) by a ik yields the solution 

(58.4) y k = a tlc y t — a ik b { . 
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Thus 

(53.5) 

dy ( _ _ „ 

d&i 11 ' 



In virtue of these equations, we see by examining (9.1) 
that the distinction between eontravariance and co¬ 
variance has disappeared* Accordingly we shall write all 
indices as subscripts on condition that we allow only 
transformations of the type (53/2) subject to (53.3)* We 
have already anticipated this in our notation y it 6 iS and 
a i3 . However we may sometimes wish to adopt a cur¬ 
vilinear coordinate system such as spherical polars. It is 
then necessary to reintroduce the distinction between 
eontravariance and covariance. This will be indicated by a 
return to the coordinates 

The transformation (53*2) is equivalent to the com¬ 
bination of the two transformations y € — y\ + b t and 
Vi — a uVi - The transformation = y't + 6* merely 
defines a translation to new parallel axes. The transfor¬ 
mation 

(58.6) y'i = 

subject to the six conditions (53.3) is said to define an 
orthogonal transformation* It follows from (53.3) that 
the determinant ] a ti \ of an orthogonal transformation 
is either + I or— 1 and we say that (53.6) defines a 
positive or a negative orthogonal transformation respec¬ 
tively. It is well known* that the y * axes form a right- 
handed or a left-handed system if the orthogonal trans¬ 
formation is positive or negative respectively. Further a 
positive orthogonal transformation defines a rotation of 
the axes about the origin. 

An alternative set of equations to (53.8) is obtained 
by considering 

* W* H* McCrea, "Analytical Geometry of Three Dimensions** 
pp* 10—13. 
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ds 2 = dy i dy i = a H a ki dy f dy k = i §k dy $ d§ kf 
from which we deduce that 
(53.7) ^a a ki = 

Ex. Show that in tiie rotation of axes defined by the positive 
orthogonal transformation (5G.G), is the cosine of the angle 
between the #1 axis an d the y f axis. The components a if are 
therefore the direction cosines of the y r t system with respect 
to the y t system* 

§ 54-. Rotations 

Equations (53.2) may be interpreted from another point 
of view. We could say that they transform the point 
P whose coordinates are y t into the point P whose coor¬ 
dinates are y t referred to the same system of rectangular 
cartesian axes* We then call (58.2) an affine transfor¬ 
mation** If in addition conditions (53.3) together with 
\a if \ = l arc imposed, we have the most general rigid 
body motion consisting of a rotation followed by a trans¬ 
lation. 

Next we wish to obtain the orthogonal transformation 
which represents a right-handed rotation through an 
angle ip about a line through the origin whose direction 

* Affine Geometry. According to Klein’s Erlangen Programm, 
(F* Klein, Math. Ann* Vol 48* (1893), p, 83), a geometry comprises 
a system of definitions and theorems that express properties which 
arc invariant with respect to a given group of transformations, For 
example, if the group of transformations consists of all rigid body 
motions, namely translations and rotations, then the geometry is 
termed metrical. This is the geometry of Euclid, less the similarity 
theorems, and its most important concepts arc distance and angle* 

Let us now examine affine geometry, defined by the group of 
transformations (53*2). Suppose A 4 to be the vector joining the point 
with coordinates to the point with coordinates y t * Then 
At = y t — z t% and we immediately deduce from (58*2) that the 
transformation law of vectors is A f = a if Aj* We now choose two 
parallel vectors A t and B t , The requisite conditions are 
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cosines are . By a right-handed rotation, we mean that 
a right-handed corkscrew will move outwards from the 
origin along the direction l t when it is rotated through 
an angle ip about this direction. Using the vector notation 
we have (fig.) 



OP = 0P + PL + LP ; 

where Q is the perpendicular from P on the line through 
the origin whose direction cosines are l i7 the angle 


AJB l =r AJB 2 = But each fraction equals + a lt 

a- AjS lt and similarly equals 
AJB* and A^fB^ Thus we have AjB l — AjB^ — AjB a . That 
is, the vectors A t and B t are parallel* Hence the parallelism of 
vectors is invariant in affine geometry. 

Further we see from (58,2) that every finite point transforms 
into a finite point. Therefore the plane at infinity is invariant, and 
so we can distinguish between a non-central and a central quadric 
according as the plane at infinity docs or does not touch the quadric. 
However we cannot define either angle or distance in affine geo¬ 
metry, because they are not invariants under the affine group of 
transformations (53.2), 
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PQP — QP = QP, the plane PQP orthogonal to 
and PL is the perpendicular from P to PQ. Since 

OP=Vu we liave OQ — ltl k y k and thus PQ = khy k —Vi - 
Therefore PL — (1—cos ip)PQ = (1 — cos ip)(lihy k — Vi)* 
Further LP is orthogonal to both PQ and l t in such a 
way that PQ , l t and LP form a right-handed system. 

Then by (39.1} the unit vector in the direction LP is 
given by — Vi)hlPQ = — emViklPQ* Hence 

LP = PQ sin ip (— e iJk y,l k IPQ) = e iik l j y k sin ip. Thus we 
have 

Vi = Vi + ( l — eos fWthy* — Vi) + sin ip e m l } y k7 
which can be written 

(54.1) y t — n ik y k 
where 

(54.2) a ik = cos ip8 ik + (1 — cos ip)l(l k + sin ipe m l jt 

In the theory of elasticity we shall be particularly 
interested in infinitesimal rotations* in which case 
cos ip I and sin ip ^ ip. The infinitesimal rotation is then 
represented by 

Vi = Vi + *pe m l t y k * 

That is* 

(54.3) Vi = y t + s tk y k , 
where 

(54.4) s ik = h • 

It is clear that s ik is anti-symmetric. Conversely, let us 
consider equations (54.3) given that s ik is anti-symmetric. 
Then equations (54.4) comprise only three equations fur 
the three unknowns l i7 whose solution is in fact — B^jip y 
l 2 = — s m lip and l 2 = — s^/ip. Therefore equations (54,3) 
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always represent an infinitesimal rotation if tp is in¬ 
finitesimal and s ik is anti-symmetrie. 

Ex. Calculate corresponding to a rotation of 00 5 about the y t 
axis. 

g 55, Cartesian tensors 

A Cartesian tensor of the M-th order in a three- 
dimensional Euclidean space is defined as a set of 3* f 
quantities which transform according to equations (9*1) 
when the coordinates undergo a positive orthogonal 
transformation. This is a less stringent condition than 
that imposed on a tensor. So we see that all tensors are 
Cartesian tensors but a Cartesian tensor is not necessarily 
a tensor in the usual sense. In virtue of (53,5) we have 
that " s a Cartesian tensor of the Af-th order 

if the transformed components satisfy 
(55.1) i = k dr k fit k A. k t t j 

on change of the coordinates by the positive orthogonal 
transformation = a iS y jw We sec from (53,6) that both 
y f and their differentials dy t arc Cartesian vectors. Also 
the Kronecker delta is a Cartesian tensor of the second 
order because 

= a i, a iT = Sit 

in virtue of (53.7), Further we deduce from (38.2) that 
the permutation symbol e m is a Cartesian tensor of the 
third order. Thus a ik and s ik introduced in the last 
section are Cartesian tensors of the second order. 

The fundamental tensor of the Euclidean space is the 
Kronecker delta <5^, Hence all the Christoffel symbols are 
zero and so the comma notation for covariant derivatives 
now denotes the familiar partial derivatives which are 
Cartesian tensors. 

The use of (55.1) instead of (9.1) shows us clearly that 
the quotient law of section 12 still applies to Cartesian 
tensors* 
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g 56. Infinitesimal strain 

Consider a body which is strained by the action of 
applied forces. The particle at the point P with coor¬ 
dinates y t referred to a rectangular cartesian system is 
displaced to the point P witli coordinates y { + u { * Similar¬ 
ly the particle which was at Q with coordinates is 
displaced to the point Q with coordinates z t + v t « We 
define the extension of the straight line joining the 
unstrained points P and Q to be the increase in length 
per unit of length due to the strain. That is 


(56.1) 

We have 

and 


e (PQ) = 


PQ-PQ 

PQ 


PQ 

PQ 


i. 


(PQ? = (y i -^ i )(y i -z i ), 


{ PQ? = (yt + «i — ~ vt)(y t + — Zf —1><) 

= (yt — zt)(yt — *() + 2 (y f — z,){Ui — «*) 

+ ( w t — *0{ M ( — l ’() 


= (PQ? { 


2l t (u t — v f ) (u t — v t )(u ( 


'* ) 


1--L-_ 

PQ + (PQ) 

where = (% — y t )jPQ arc the direction cosines of the 
unstrained line PQ. Hence 


e {PQ) — 


2/ r (u < — v t ) (Ui — vAiUi — Vi)]^ , 

PQ + {PQ)* I ' 


We confine our attention to the case of infinitesimal 
strain* where the assumption is made that the displace¬ 
ment vectors u t and are small compared to PQ . 
Neglecting quantities of higher order than the first, we 
obtain 


€ ( PQ) = 


h(Uj — Of) 

PQ 


Let us now choose Q to be in the neighbourhood of 
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P so that y t — z ( is small. Then by Taylor’s theorem for 
a function of three variables* 


v { = «, + (z t — y } )Uij + terms of higher order in as,— y f . 

Neglecting terms of higher order than the first in as,— y it 
we obtain that the extension e at P in the direction 
determined by the unit vector l t is given by 

lj{Zj — Vs) 


e = 


PQ 


u i t l — * 


Introduce the symmetric 
e iS by the equations 


Cartesian strain tensor 


(56,2) e i} — £( U <J + u i t i\ 

(there is no confusion between e 4i and the permutation 
symbols of the previous chapter which are distin¬ 
guished by Greek indices), and we finally obtain the 
extension e as the quadratic form 

(56*3) 6 = *?( il'ilj' 


The dilatation or expansion 6 is defined as the 
increase in volume per unit of volume* That is 


(50*4) 


AV^AV _AV 
AV AV 


where AV denotes the strained volume corresponding to 
the volume AV '■ But 


AV = d(y £ + Ui) 
AV d{ yi ) 


1 + U l,l + »2.B + %,3 + termS ° f 

higher order* 


Thus in virtue of (50*2) and (56*4) we have 
(56*5) 8 =e it . 

The components of the strain tensor are not entirely 
arbitrary* To prove this, differentiate (56*2) twice and 
obtain 


* R* P* Gillespie, Partial Differentiation, p, 60* 
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€u,ki = 

from which it immediately follows that 

( 56 * 0 ) €ii t kt + e kt t ii — e ikjt — ~ {J - 

There appear to be 81 of these compatibili ty equations* 
In actual fact some are repeated due to the symmetry of 
the strain tensor and others are satisfied identically* The 
reader is asked to verify that only six of these equations 
are independent. 

To conclude this section, we investigate several im¬ 
portant examples of strain. 

(1) Uniform dilatation* Consider the displacement 
vector u 4 — cy £ where e is a constant* We have e if — c8 ti 
and the dilatation 0 = 3c* Thus the extension at any 
point in any direction is constant and equal to one-third 
of the dilatation* 

{2} Simple extension* Consider the displacement 
vector Ui = cliljy j where c is a constant and l t a unit 
vector* We have Ufj — cl^ and so e iS — cl 4 l 4 and the 
dilatation 0 = cl ^ = c* Also = <? and thus there is 
an extension at any point in the direction l ( of amount 
equal to the dilatation* The extension in any direction 
orthogonal to l { is easily seen to be zero* If c is negative 
we refer to it as simple contraction* 

(3) Shearing strain* Consider the displacement 
vector u i — 2cl i ?njyf where c is constant and both l t 
and m { are unit vectors* A brief calculation yields 
e i} = + i*m ( ) and the dilatation 0 = 2d i m it Conse¬ 

quently the dilatation is zero if the directions l { and m £ 
are orthogonal* 

Ex* Show that a simple extension along any direction together with 
an equal simple contraction along an orthogonal direction is 
equivalent to a shearing strain along a direction bisecting the 
angle between the given directions* 
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§ 57. Stress 

The forces acting on n body are either external or 
internal. The external forces may consist either of body 
forces such as gravity which act on every particle of it, 
or of surface forces which act on the external surface of 
the body, for example the pressure between two bodies 
in contact. If j V i denotes the body force vector per unit 
volume, then the force acting on an clement of volume 
AV is FjAV* Similarly if Tj denotes the surface force 
vector per unit area, then the force acting on an element 
of surface AS is T t AS, In order to discuss the internal 
forces, we select a small element of area AS inside the 
body and denote the direction cosines of the normal to 
this element, which is approximately planar, by %. We 
call one side of the element AS positive and the other 
side negative. Then the action of the positive side on the 
negative side is the internal surface force T } AS where 
T j is the force per unit area on the element AS. It is 
called the stress vector and is in general a function of 
the coordinates of the point which determines the position 
of the element AS and of the direction cosines n s of the 
normal to AS, It is well to emphasise that T i is not 
necessarily co-directional with %. At all points on the 
external surface of the body T j becomes the external 
surface force. 

Consider a small rectangular parallelepiped with vertex 
at the point P whose edges are parallel to the coordinate 
axes. Wc form three stress vectors T {2 )j and T^ 

corresponding to the small elements of area through P 
which are parallel to the coordinate planes. The stress 
vector Ta )j will be called positive if it acts in the positive 
direction of the y t axis provided that the external normal 
is co-directional with the positive y t axis. If, however, 
the external normal is co-directional with the negative 
y t axis, then the stress vector is positive if it acts 
in the direction of the negative y f axis. In other words, 
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a stress which tends to stretch will be regarded as positive, 
whilst a stress which tends to compress is regarded as 
negative. We define nine quantities E i} by the equations 

(57.1) E ti = T m 

and we shall show that E if is a Cartesian tensor, called 
the stress tensor. 

Construct the small tetrahedron PA t A 2 A 3 such that 
the edges PA t are parallel to the y t axes. The forces 
acting on the tetrahedron are tfie body force F f AV> the 
surface forces E^ASi {no summation over i) on the 
faces opposite to and the surface force TjAS acting 
on the face A t A 2 A %l where AV is the volume of the 
tetrahedron, AS t is the area of the face opposite to A { 
and AS is the area of the face A x A t A z . Let the positive 
direction of be that of the normal drawn outwards 
from the tetrahedron and whose direction cosines are % 
and let p be the perpendicular distance from P to the 
face A y A % A z . Then AV = %pAS and AS t — n f AS. The 
equations of equilibrium obtained by resolving forces 
parallel to the axes are now 

F i AV—E ii AS i + T s AS = 0, 

where the E if occurs with a negative sign, because the 
external normals arc in the directions of the negative 
axes. We now substitute for AV and AS if then divide by 
AS and proceed to the limit as the tetrahedron shrinks 
to zero, in which case p tends to zero. The result is 

(57.2) Ti = 

It follows from the quotient law that E t} is a Cartesian 
tensor, and we can by (57.2) calculate the stress vector 
at any point, corresponding to any direction at that point 
in terms of its direction cosines and the stress tensor. 

We now cite several important eases of stress 

(1) Normal Stress. The vector Ti is co-direetlonal 
with Az n i ■ We see from (57,2) that E if — C6a where C 
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is a constant. Hydrostatic pressure is an example of 
normal stress for which C is negative. 

(2) Simple Tension. Consider the stress tensor 
E iS — Clilj where C is a constant and l t is a unit vector. 
Then the stress vector in the direction i* is T t —Cl t 
and it is thus co-directional with ± However the stress 
vector in the direction orthogonal to l € is T i — ^ro* — 0. 
If C is negative the stress is called a Simple Com¬ 
pression, 

(3) Shearing Stress, This is specified by the stress 

tensor E H — where C is a constant and 

l t and nii are unit vectors. 

Ex, Show that a simple tension along any direction together 
with an equal simple compression along an orthogonal direc¬ 
tion is equivalent to a shearing stress along a direction bisecting 
the angle between the given directions. 

§ 58. Equations of equilibrium 

Let us consider a body of volume V in equilibrium, 
which is enclosed by the surface S , We resolve forces 
parallel to the axes and obtain 

j y + J g T J dS = 0. 

Substitution from (57.2) yields 

j y F j dV -j- | s E if n,dS = 0 , 

which on application of Gauss’s theorem * becomes 
j v F,dV+ j y E„' i dV=0. 

That Is, 

j v (F, + E Uii )dV= 0. 

This equation is an identity, being true for any volume V , 
Hence 

(58.1) F, + E iKi = 0. 

* D, E, Rutherford, Vector Methods, p. 74. 


§ 59 CARTESIAN TENSORS - ELASTICITY 95 

The moments of a force about the axes are the com¬ 
ponents of the vector product of the force vector and the 
position vector of any point on the line of action of the 
force. Thus in tensor notation, the moments of the force 
F| about the axes arc represented by the vector e iik y t F k . 
Now take moments about the axes for our body in 
equilibrium and the result is 

E k dV + ^ s e m y $ T k dS = 0 . 

We substitute from (57.2), apply Gauss’s theorem and 
obtain 

j v F k dV + j v {e iik y f 

That is, 

j v e ifiVAFk + E lkl )dV + j y e m d n E lk dV = 0 . 

The first integral is zero in virtue of (58.1) and we have 

jy e itkEltdV — 0 , 

This equation is also an identity and consequently the 
integrand tiu k E lk vanishes, from which we deduce that 

(58.2) E if — E Si . 

Thus the stress tensor is symmetric and the elastic 
equilibrium equations for a body are given by (58,1). 

§ 59, Generalised Hooke’s law 

In the elementary theory of elasticity, Hooke’s Saw 
states that the tension of an elastic string is proportional 
to the extension. In other words, stress is proportional 
to strain. The corresponding assumption in the general 
theory of elasticity is that the stress tensor is a linear 
homogeneous function of the strain tensor. That is 

(59.1) E u — c im e kl . 
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It follows from the quotient law that c im is a Cartesian 
tensor of the fourth order, and it is called the elasticity 
tensor* Further from the symmetry of E u and e kt we 
find that c im is symmetric not only with respect to the 
indices i and j but also with respect to k and 
A body is said to be homogeneous if the elastic pro¬ 
perties of the body are independent of the point under 
consideration. This means that the components of the 
elasticity tensor are all constants for a homogeneous body* 
We call a body isotropic if the elastic properties at a 
point are the same in all directions at that point. This 
means that the elasticity tensor c iikl transforms to c im 
itself under any rotation of axes. 

| 60. Isotropic tensors 

A Cartesian tensor which transforms into itself under a 
rotation of axes is called an isotropic tensor. We have 
already met two isotropic tensors, namely 6 H and c iik . 
We shall now search for the most general isotropic 
tensor c im of the fourth order. Its transformation law 

(55.1) becomes 

( 60 . 1 ) C ijkl = fti r ®j a a k t a iU C TBtu* 

Let us rotate the axes through 180° about the y z axis. 
We deduce from (5-4.2) that a ik =— + But 

for this transformation ^ — 0 and 1$ = 1, and so the 

only non-zero components of a ik are 

#11 = —* 1 , ” 1 ) #33 mm Hr L 

By direct substitution in (60.1) we obtain c im — — c iiki , 
that is, c im — 0 in the following cases:- 

(1) any three of the indices equal to 1 and the other 
equal to 8. 

(2) any three of the indices equal to 2 and the other 
equal to 8. 

(8) any two of the indices equal to 1, another equal to 
2 and the other equal to 8. 
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(4) any two of the indices equal to 2, another equal 
to I and the other equal to 3. 

Similar results are obtained by the corresponding 
rotations through 180° about the y 1 and the y 2 axes. 
Hence the only components which survive are those 
where the four indices arc equal or equal in pairs. 

Let us now rotate the axes through 00° about the y% 
axis. For such a transformation w p e deduce from (54.2) 
that a ik = l t l k + e m l s and so the only non-zero com¬ 
ponents of a ik are 

#12 = — 1; #2! = + I; # 33 « + L 
By direct substitution in (60.1) we discover that 
c nn = cam* 

c naa — c mit — £2233* c 33 ii = £3322, 

Cl212 = ^2121 1 C1313 — C £323 , C3131 — e 3 232 * 

C1221 = Cana , cim = £2332* 03113 — C3223, 

Similar results are obtained by the corresponding rotations 
through 90° about the y t and the y 2 axes. Hence we can 
collect all our results in the form 

c iiii ~ 

c tiii — c iikk — c llii ~ 
c Uii = c ikik == c lHi = 
c ifii — Cikki “ c tm = C lkktt 

where i t j> k and l are unequal and the summation con¬ 
vention does not apply. All other components are zero. 
The most general solution of equations (60,2) is then 

(60.8) c im — Xd (j d kl + pd ik d f i + + K &mt* 

where A, pt, v and k are Cartesian invariants, and S im = 1 
if all four indices are equal and otherwise zero. 

Finally we carry out a small rotation, which is re- 


(60.2) j 


H 
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presented by a ik = d ik + $ ik9 where $ ik is skew-symmetric 
and of the first order in small quantities. Substituting 
in (60*1) and retaining only terms of the first order 
we have 

9 ir C rm + S i* C i*kt H“ S M C im + S tu € iiku = 

Select i — 2, j — k = l — 1 . Then in virtue of s 12 — —s 2l f 
the non-vanishing terms of this equation yield 

c un “ C 22U + c 312t + C 2112 ■ 

Substitute in this from (60*8) and the result is 

A + ^ + v + A + /* + ^ 

which gives k = 0 , so we can now write 
(60.4) c tm — k&ij&ta + ft&ik&H + vditdit* 

This tensor is obviously isotropic and so it is the most 
general isotropic tensor of the fourth order* 

Ex. 1. Prove that the most general Isotropic tensors of the second 
and third orders are M iS and A%* respectively where A is 
an invariant. 

Ex. 2* Prove that e ml defined by (G0.3) satisfies the symmetry 
relations and 55 * 


§ 61. Homogeneous and isotropic body 

In this section we confine our attention to a body 
which is both homogeneous and isotropic* Then the 
elasticity tensor is isotropic with constant components. 
Wc therefore require an isotropic tensor c mi which is 
symmetric in both i and j and in k and i, and whose 
components arc all constants. To obtain this tensor, sub¬ 
stitute from (60.4) into the equation c im = c mk and we 
obtain 

(// — * 0 ( 5 ** — &u 5 ^) = 0 . 

In this equation put i = fc = 1,7 = I = 2 and we see that 
jM = v* Thus the isotropic tensor 

(61.1) c mi — ^5 w 5* ( + ft(6 ik d n + 5 n 5 Jt ) 
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satisfies the required symmetry. The generalised Hooke’s 
law (59.1) for a homogeneous, isotropic body now 
becomes 

Ea = [A5**5^ +/i(5*jfc5*| + 8 n 6 ik )]e kl 
where A and are constants . This equation simplifies to 

(61.2) E ts — A86 ft + 2 

where 0 is the dilatation defined by (56.5). Contracting 
i and we obtain the equation 


(61.8) 0 = Ea = (BA + 2fi)0 


which connects the two invariants 0 and 0 , 

We now solve equations (61.2) for the strain tensor 
in terms of the stress tensor and in virtue of (61.8) we 
obtain 




It is usual to work in terms of Young’s modulus E and 
Poisson’s ratio a defined by 


(61.5) E - ^ (8A + 2fl) a 

A + n 

Equation (61.4) then becomes 


A 

2(A + fi) 


(61.6) Cft — — { a&6 ft + (1 + 0 )Eft}, 


We can obtain stress compatibility equations by sub¬ 
stituting this expression for e iS in (56.6). 

Sometimes we require the equations of equilibrium in 
terms of the displacement vector u *. To obtain these we 
substitute from (61.2) in (58.1) and obtain 

F f + Ad 4 dft + — 0 

which becomes by (56.2) 

Ff + A0 f + + u i a) — 
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But 0 — % = i(u ti + u ifi ) = u ifi * Consequently 
Oj = u t ii9 and the atiove equations take the form 

(6L7) F i + (X + p)Q wi + pV*u § = 0 

where V 2 is the Laplacian operator. Now that we have 
reverted to displacements u t we do not require any com¬ 
patibility equations. 

Ex. Show tluit if an Isotropic, homogeneous body is in equilibrium 
under no forces, then — 0. 

§ 02 . Curvilinear coordinates 
Many problems in elasticity can be investigated more 
conveniently by means of curvilinear coordinates x* in 
which the line-element ds is given by d$* = g ii dx i dxK In 
this case we must return to the distinction between 
contravarianee and covariance. 

The coordinate system serves only to describe the 
actual strains and stresses. The laws of elasticity are 
themselves independent of the coordinates. Therefore 
these laws can be formulated as tensor equations. We 
recall (section 9) that if a tensor is zero in any one co¬ 
ordinate system, it is zero in all coordinate systems. 
Consequently if we write down tensor equations, which 
reduce in cartesian coordinates to the results already 
established, then these tensor equations express the 
theory with respect to any curvilinear system. Ac¬ 
cordingly we immediately verify that the laws of elasticity 
are represented by the following tensor equations:- 

e = e ti l*V, e ti = Hu tif + u f i ), 0 = g {i e„, 

e ii,kl + ^kl,H - e ik,il - e H,ik — 0, 

T t = E, f n>, E tf = E h , F* + = 0, 

Eif = c ( yitje* 1 , c iikl = c fikl = c iilk *. 

* It can also be proved that c kUi — c im , See I. S. Sokolnikoff, 
Mathematical Theory of Elasticity, section 26. 
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If the body is homogeneous and isotropic, we have in 
addition 

e im = tyiight + MSkigif + 

E iS = X0g £i + 2*pe U9 0 = g** Eh — (3A + 2/*)0, 

e a = jg {— v@gij + t 1 + 

F s + (X + p)d tf + (ig™u Ura = 0 . 

In all these equations I 1 is a unit vector specifying the 
direction of the extension e, n* is the unit vector normal 
to the small element AS, and commas once more denote 
eovariant differentiation. 

It is important to note that the components of the 
vectors u it T t and F t may not possess a physical signi¬ 
ficance of the correct dimensionality. For example, the 
exercise of section 5 shows us that the second component 
of the acceleration vector in polar coordinates is an 
angular acceleration. It will suffice to discuss the force 
vector F& whose components in a cartesian system 

are F £ = F* (say). Then F t = ~ F £ and F* — ~ F* 

cat* 0$ 

are the covariant and contravariant components in the 
x i system, where we have put x 1 — y £ the cartesian 
variables. The component of F\ which is the physical 
force vector, at any point in the direction of the 
unit vector I 1 = l £ is the invariant F i l i = F i l i = 
g ij F i lK The contravariant unit vector in the direction 
of the x 1 axis is dJ/VIiI* Therefore the actual physical 
components of force along the x l coordinate curve arc 
gii^/VST* In *h e exercise of section 5, 
we have g u — 1 , g 12 = 0 , g 22 — r 2 and so the physical 
components of acceleration in polar coordinates are 
d*r (dd \ 2 , d*Q t dr dO .. . 

W ~ r Uj antl r dF + 2 JtTt’ whlch rcprcsent the 

radial and transverse accelerations respectively. 
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Similarly the tensors e i} and have no direct physical 
significance. Consider the strain tensor e i} which is 
connected with the cartesian strain tensor bv 
dx k dx l 

€ u = d9 * kl ' The P h y* ical com P°nent of the strain 

tensor associated with the directions of the unit vectors 
P and m * at a point may be defined to be the invariant 
PmK 

As an example, let us discuss cylindrical coordinates 
r, 0 and z given by 

Vi = r cos 0, y 2 — r sin 0, y 3 = z. 

The non-zero components of the fundamental tensor arc 
gu^ 1 - = g» = l 

from which we deduce that 


g 11 _ i f g« = i, g« = 1 ( gi/ = 0 if i ^ j. 

So the only surviving Christoffel symbols of the second 
kind are 



The actual physical components (u a} u fi , u 7 ) of the dis¬ 
placement vector u t are 

1 

U a Uj ) UjJ “ Uy — l&j. 


Further calculation will show that the components of 
the strain tensor are 
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whilst the physical components associated with the direc¬ 
tions of the coordinate curves are 


du a 


1 dufl . u a 
r 


du y 


"IF' e "~TW + T’ 6yy ~~W’ 


(1 Bu a duf uA _ , (dtiy duA 

e * fi ~ * \T W + 1 F~T) ’ gay ~ 4 [IF + IF)' 
e . _ i (ta + I tl\ 

eyf ~ i \dz + r 90)‘ 


The dilatation is 


o „ 3w» , 1 , Suy , u a 

e = t“«- w + 7W ' + ^ + T' 


Ex, Find the physical components associated with the directions 
of the coordinate curves of the strain tensor in terms of the 
physical components of the displacement vector when the 
coordinates are spherical polars* 


| 68. Mechanics of continuous matter 
Let us discuss the motion of a continuous medium by 
the Eulerian method. Instead of following the path traced 
out by a particular particle, we focus our attention on a 
definite point P of the medium whose coordinates referred 
to a cartesian coordinate system are y { . Let us denote 
by Uf the velocity vector of that particle which happens 
to be at P at time t Then u ( is a function of and t 
After a further interval of time At has elapsed, the particle 
which was originally at P is now at the point y t + u t At 
with velocity u t + Au { * Hence u { + Au t is the function 
u ( at the point y t + u { At and at time t + At On expand¬ 
ing by Taylor’s theorem we have 

0 a 

U 4 + Aui = Ui + UjAitiij 4- At * 

The acceleration vector / f at P is the limit of AuJAt as 
At tends to zero. We therefore deduce that the acceleration 
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vector is given by 

( 68 . 1 ) /, = u,u iti + ^. 


§ 63 


Now consider a volume V of the medium, bounded by 
the surfaced. The massitf contained inside V\sM=$ y pdV, 
where the density p is a function of y t and L Thus the 

rate of increase of mass is ~ = J ^ dV. Let n £ denote 

the direction cosines of the external normal to the small 
element of surface AS . Then the rate of mass flow out¬ 
wards across that element AS is pr^UiAS. Hence the rate 
of increase of mass is also expressed by the surface 
integral —J sp^^dS. Consequently we have 

pngUidS + J* ~ dV =* Q. 

We apply Gauss’s theorem to the surface integral and the 
result is 

f^l'dv+j^ir- o. 


This equation is an identity and so we derive the equa¬ 
tion of continuity 


( 68 . 2 ) (/**,)., + = o. 


Physically this equation expresses the principle of con¬ 
servation of mass. 

In section 58 we discussed the equilibrium of a con¬ 
tinuous medium. The equations of motion can be deter¬ 
mined in the same way. If we replace F i by F s ~pf $ 
these equations become 

(63.8) pft=F i + E iiti . 

We now substitute for f, from (63.1) and obtain 
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pUfUf i + p = F f + E if i , 

This equation can be written, in virtue of (68.2), in the 
form 

9 

(63.4) QW|«i — E ti ] iX + ^ {puj) = F ,. 

(63.2) and (63.4) constitute the equations of motion of a 
continuous medium. 

When the coordinate system is curvilinear the equations 
(03.2) and (63.4) may be expressed in the tensor forms 

(08.5) (pu‘) >( + ^ = 0, 

(63.6) (pu*v? - W*) ti + | (put) = F*. 


Solutions 


p, 88. Ex. The only non-vanishing components of are 

^ii “ '— h ®n = + lj — I - 

_ dtl<* l bus tin 

p. 108. Ex. e aa = _ ; m = T M + T S 


fttt-L 


C YY 


Ua COt 0 

. + — +- n/n 

rsintJ btp r r 


1 / 1 bn# 1 bUy COt0 \ 

efiy ~ t \ 7 toe v + v w ~ ~r M T 

2 \rsin 0 btp br r f 

S'). 

2 (r M br r ) 













CHAPTER VIII 


THEORY OF RELATIVITY 

§ 64. Special theory 

In classical mechanics, the position of a point in space 
at which an event occurs can be determined by its three 
space coordinates x\ x®, x 3 referred to some rectangular 
cartesian system. Also an observer can measure the time t 
at which the event takes place by means of a clock. An 
event is then fixed in both space and time by the system 
S which is comprised of the four numbers x\ x\ x 3 and L 

Einstein examined the concept "simultaneity* and came 
to the conclusion that "simultaneous events at different 
points* has no meaning without further qualification. 
Continuing his study of fundamental ideas, Einstein 
arrived at the Special Theory of Relativity which he 
based on the two principles: (1) it is impossible to detect 
the unaccelcraied tramlutory motion of a system through 
space t ( 2 ) the velocity c of a ray of light is a constant which 
does not depend on the relative velocity of its source and 
the observer . 

Let us now consider two systems S and 8 which coin¬ 
cide at the time t — 0 , such that S moves with constant 
velocity V along the x 1 axis of the S system. Then the 
Lorentz Transformation, which can be deduced from 
the two principles of special relativity, connects the 
space coordinates and the time of both systems by the 
equations 

x 1 — l Jfx 1 — Ffc), = x® x 3 s= x 3 , i = {}{t — Vx l fc 2 ) t 

where /? — (1 — F 3 /c®)~^ ( We easily verify that 
— (dx 1 )* — (dx 2 )® — (dx 3 )® + c*{dlf 
= — (dx 1 ) 2 — (dx®J a — (dx 3 ) 3 + c®(d*j* 
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The invariance of this equation with respect to Lorentz 
transformations suggests that the Minkowski space 
defined by the metric 

(64.1) da 3 = — (dx 1 ) 2 — (dx®)® — (dx 3 )® + c*(dx *) 3 

where we have written x 4 = f, is appropriate for the 
geometrical discussion of special relativity. We denote the 
line-element of this four-dimensional space by da {not 
by ds) in order to emphasise that da is not the physical 
distance between two neighbouring points. 

The Minkowski space is flat and its signature, which 
equals the excess of the number of positive terms over 
the number of negative terms in its metric, is — 2 . It 
is well-known* that there is no real transformation of 
coordinates which will reduce (64.1) to the metric of a 
four-dimensional Euclidean space, whose signature is 4. 
Thus the geometry of Minkowski space differs in many 
respects from Euclidean geometry; for example, there 
exist real null curves (see (15.2)) and real null-geodesics. 

In this chapter Latin indices will have the range I to 3 
whilst Greek indices will range from 1 to 4. The velocity 
u of a particle, which is at the point x i has the components 
u* = dm*jdt referred to the system S . It follows from 
(64.1) that 



The four-dimensional Minkowski momentum vector 
dx a 

is defined by m Q c -where m 0 is a constant. The special 
da fjfa* 

theory identifies the fourth component m 0 c— with the 

mass m of the moving particle. In virtue of (64.2) we 
have 

/ u 2 

(64.8) m —m 0 ll —^ 

* W. L, Ferrar, Algebra, p* 154. 
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The constant m 0 is the mass when u = 0 and so it is 
called the rest-mass of the particle. The mass m, which 
clearly increases with the velocity, is called the rela¬ 
tivistic mass of the particle. The components 


dx* dx l di dx i 

m ° c ~dtdB = m ~dT 


and are evident generalizations of the Newtonian mo¬ 
mentum vector. 

We define the four-dimensional Minkowski force 
vector F* by 


(64,4) 


d I dx*\ 

= di{ m lu) 


The Newtonian force vector is X* ■ 


and so 


We obtain by expansion from (64,8) that 


me 2 = m 0 c 2 + 


1 

2 


m 0 u 2 + *,, 


and so the special theory identifies the energy E as¬ 
sociated with a particle by the equation E = me 2 . 
Therefore 

l 1_ ?J = t- 

The motion of a particle which moves under the action 
of some force system can be represented in Minkowski 
space by a curve, called the world-line of the particle. 
If no forces act on the particle, we see from (64,4) that 
dPxpjda 2 — 0, Thus the world-line of a free particle is a 
geodesic of the Minkowski space. 

The velocity of a light ray is the constant c, and so 
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we see from (64.1) that for such a ray da — 0* Accor¬ 
dingly the world-line of a light ray is a null-geodesic of 
the Minkowski space. 

In order to discuss the mechanics of a continuous 
medium, we introduce the symmetrical four-dimensional 
energy-momentum tensor 2 W defined by 

T u = T fi = pu l u* — E**; T a = T Ai = pu { ; - p, 

where p is the density and E (i is the Cartesian stress 
tensor defined in section 57, Then the special theory 
generalises (63,5) and (63.6), which are the equations of 
motion of a continuous medium into 

(64.5) Tg = F*. 

If we change to spherical polar coordinates r, 0 and ipj 
the metric of Minkowski space becomes 

(64.6) da 2 = — dr 2 — r 2 d6 2 — r 2 sin 2 Odip 2 -j- c 2 dl 2 . 

§ 65, Maxwell's Equations 

The classical theory of electrodynamics*, according 
to Lorcntz, is specified by the electric potential <p which 
is a scalar and the magnetic potential A £ which is a 
vector. The electric field strength vector E { and the 
magnetic field strength vector Il t are derived from 
these potentials by the equations 

r. , 1 dAi 

E i = -g rad*—*-. 


H { — curl A f , 


* M, Abraham—R* Becker, Electricity and Magnetism. 

C, A, Coulson, Waves, chap, VII. 

D. E, Rutherford, Vector Methods, p. 126. 
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Using electrostatic units, Maxwell’s equations are 


(65.1) 


div Ei = 4jrp, 
div /I* = 0 , 


i * f 1 dH t 

curl e ‘ + 7 St 

1 IT 1 $ E i 

curi// '-7 3 T 



where /, is the current density vector and p is the 
charge density. 

In Minkowski space* with the metric (64.1), let us form 
the four-dimensional potential vector 0 a and the four- 
dimensional current density vector J a defined respec¬ 
tively by 

0a ■= ( * -^2 i ^ 3 1 mp*)* 

Jet = (/is /as ?3> p)s 

with respect to a particular coordinate system. Next we 
introduce the skew-symmetrical tensor defined by 


V&0 — *&0 r aL 


d0* d00 

dx$ dx* 


and w'e immediately calculate that its non-vanishing 
components in the given coordinate system are 

^23 = ^?32 ” H i; — ' ^21 — 

Vl4 ” Vn = cE l* *?24 “ ^?42 = cE 2> *?34 == ^43 = 

The non-vanishing contravariant components may 
now be obtained and arc 


= H x ; = Th; 71 *=—*}* = H z ; 

v u = - n a =— =—}?«=— T) M =~r) a =— 

c c 0 

We now write Maxwell’s equations (65.1) in terms of 
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rj and J and the results are readily verified to be res¬ 
pectively 

dif 2 di j 43 in 

dx 1 dx 2 dx 3 c 9 

d??23 , , 57?i3 „ 

3a : 1 " l " 3a:* ^ 3a? ’ 

\ d M ' dr!ii - » 

3a * ' dx* i " dx* ~ * 
drf 1 dyf* drj i3 drj H _ 4 tt { 

dx 1 dx 2 dx 3 dx 4 c 

The first and last of these equations combine together 
into the form 

(65.2) ^ J«, 

whilst the remaining two are accounted for by the 
equations of the set 

( 6 5,8) " 1 “ Vfly** + — 0 

which do not vanish identically. 

We have accordingly written Maxwell’s equations in 
tensor form in Minkowski space. Thus they are invariant 
under the Loren tz group of transformations. 

§ 66. General theory 

We now turn to the General Theory of Relativity which 
was developed by Einstein in order to discuss gravitation. 
He postulated the principle of covariance, which 
asserts that the laws of physics must be independent of 
the space-time coordinates. This swept away the privi¬ 
leged role of the Lorentz transformation. As a result 
Minkowski space was replaced by the Riemannian V A with 
the general metric 

(66.1) da 2 — gafidxPdxP* 

Einstein also introduced the principle of equivalence, 
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which in essence states that the fundamental tensor 
g a ft can be chosen to account for the presence of a gravi¬ 
tational field. That is, g a $ depends on the distribution of 
matter and energy in physical space. 

Matter and energy can be specified by the energy- 
momentum tensor which in the special theory satisfies 
the equation T*£ — FP. The only forces, namely those 
due to gravitation, are however already taken into 
account by the choice of the fundamental tensor g a ^. 
We therefore ignore FP and, in accordance with the prin¬ 
ciple of covariance, the energy-momentum tensor must 
now satisfy the equation — 0. Wc shall write this 
equation in the equivalent form = 0 where 

T* = gp y T*y is the mixed energy-momentum tensor. 
The problem now is to determine as a function of 
the g a fi and their derivatives up to the second order, 
bearing in mind that T* fia = 0. We recall from (34.3) 
that Einstein’s tensor defined by 

( 66 . 2 ) — iRS$ 

satisfies the equation G a ^ a — 0. The equations of motion 
require T a fi a = 0, but very remarkably GJ^ ft — 0 is 
an identity in Riemannian geometry. This led Einstein to 
propose the relation 

(66.3) kT* 0 + = 0. 

In effect these equations form the link between the 
physical energy-momentum tensor T a ^ and the geome¬ 
trical tensor Gf 0 of the V 4 of general relativity. In order 
that Newton's theory of gravitation can be deduced as 
a first approximation from Einstein’s theory, it was 
found necessary to choose k = STzk/c* where k is the 
usual gravitational constant 6.664 x 10 -8 cm 3 . gm.“* 
sec. -2 . The value of c is 2.99796 x IG 10 cm, sec -1 , and 
so k is 2.073 x 10“ 48 cm. -1 gm," 1 sec. 2 in e.g.s. units. 
In the special theory, the world-lines of free particles 
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and of light rays are respectively the geodesics and the 
null-geodesics of Minkowski space. The principle of equi¬ 
valence demands that all particles be regarded as free 
particles when gravitation is the only force under con¬ 
sideration. Then it follows from the principle of co¬ 
variance that the world-line of a particle under the action 
of gravitational forces is a geodesic of the with the 
metric ( 66 . 1 ), Similarly the world-line of a light ray is a 
nuil-geodesic. 


§ 67. Spherically symmetrical metric 

Genera! relativity discusses several important problems 
in which the coordinate system r, 0 , ip and t is such that 
the metric takes the form 

(67.1) da 2 — — e k dr 2 — r 2 d0 2 — r 2 sin 2 0 dip 2 + c 2 e v dt 2 > 

where X and v are functions of r. A metric of this type 
is said to be spherically symmetrical. It is a generalisa¬ 
tion of the special relativity metric (64.6), which is ex¬ 
pressed in spherical polars. The coefficients of dr 2 and di 2 
have been selected as exponentials in order to ensure 
that the signature of do 2 is — 2 , Let us write x l = r, 
x 2 = 0, ft 3 = y) and = cL Then the non-zero components 
of the fundamental tensor are 


Sn ~ ~ g 23 = — r2 > 633 = — r 2 sin 3 0, g 41 = e r . 

The determinant g becomes 


g — — r 4 sin 2 0 


and hence the non-zero components of the conjugate 
symmetric fundamental tensor are 


rll_-g22 = — 


r- 


. Jl **33 
* 6 


r 2 sin 2 0 


* g u = 


A brief calculation shows that the only non-vanishing 
Christoff el symbols of the second kind are 


I 
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(67.2) 


{AH'- {A)-f (AH' 
{A }— nrK (A)—{As} =oote - 


(A)— r * bi ’ 9 *"*• 


{14} “ *-'• 


where a dash denotes differentiation with respect to r. We 
now evaluate the components of the Ricci tensor by 
means of (33.2) and the results are 


(67.8) 


R 


n 


- -X'-Wv' + W‘ + lv'\ 


R iZ = cosec 8 61 = — 1 + e~*{l — \ r V + £/»>'}, 

R 


■M 




t R a p = 0 for a ^ ft. 

A further calculation gives us the curvature invariant 

(67.4) + 

On substituting this in ( 66 . 2 ) we obtain the components 
of the Einstein tensor for the spherically symmetrical 
metric (67.1) 

^ +r"+^+i-i • 


(67.5) 


k — 0 for a ^ /?* 


Ex. Find the necessary and sufficient conditions that a space with 
a spherically symmetrical metric be an Einstein space. 
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| 68* Schwarzschild metric 

We now seek the spherically symmetrical metric 
(67*1) consistent with the existence of one gravitating 
point particle situated at the origin, and surrounded by 
empty space* When the origin itself is excluded from our 
discussion, the energy-momentum tensor T is zero at 
all points* It follows from ( 66 , 8 ) that G— 0 , which 
yields in virtue of (67*5) the equations 



The solution of ( 68 . 8 ) is readily found to be e *= 1—3Bm/(c®r), 
where the constant of integration m introduced in this 
way can be identified physically with the rest mass of the 
gravitating particle. Further, wc subtract ( 68 * 8 ) from 
(68.1) and obtain e“ A (A' + v')}r = 0 . That is l* + v' = 0 , 
from which we have X + v = k t where & is a constant. 
Thus e v = e k {\ — 2 m/(c 2 r)}* We can now verify that 
equation ( 68 . 2 ) is satisfied identically. However, at large 
distances from the gravitating point, the metric should 
approximate to the metric (64*6) of special relativity* 
Therefore, we must select k =* 0 , We thus obtain the 
Schwarzschild metric 

( 2m\ —1 

1 — ^ I dr 2 — r 2 dO 2 — r 3 sin 3 Odip 2 

+ c! ( i -S)‘“' 

Ex* Show that a space with SchwaizschtId ’3 metric is an Einstein 
space, but not a space of constant curvature. 
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§ 69. Planetary motion 

Let us investigate the motion of a planet in the gravita¬ 
tional field of the sun. The sun will be selected as a 
gravitating particle and the planet as a free particle 
whose mass is so small that it does not affect the metric, 
and whose world-line is then a geodesic in the V 4 with the 
Schwarzschild metric (6S.4), The geodesics are deter¬ 
mined by the four equations (26,4) in which we now, 
of course, replace s by a. We shall omit one of these 
equations, in practice the most formidable one involving 
and replace it by (27.1) which is of the first order 
and is satisfied along the geodesics. When this is done, 
we no longer require the Christoffel symbols of the type 

| The remaining non-vanishing symbols of the second 

kind can be calculated from (67,2) and may be found 
to be 


1 2 1 1 

(3\ 1 

[4 \ _ m / 

2 m\ 

( 12 ) r ' 

K 

II 

r—1 

(I4j e*r* \ 

1 c*r) 


{ 3 2 8 }=- s in0 c ° S 0, cot0. 


Hence 

(69.1) 


the four equations of the geodesics arc 
d 2 0 2 dr dQ , fdw\ 2 

^ + 75;5- s ‘ n(>cos8 (J)= 0 ' 


drw 2 dr dm „ „ _ dtp d6 

do* + Tfolk + '~ t:0 °dada 



THEORY OF RELATIVITY 


117 


§ 69 


We may assume that the planet moves initially in the 
plane 0 — ?r/2. That is, dOjda and cos 0 are both initially 
zero. Then (69.1) tells us that dPdjdo* is also zero. Repeated 
differentiation of this equation shows that dW/da *vanishes 
at t — 0 for all L Hence 0 = jt/2 permanently, and the 
above equations simplify to 


(69.2) 

2 dr dip 
da 2 r da da 

(69.3) 

d 2 t 2m / dt _ 

do* c*r \ c*r) da da ~ 

(69.4) -| 



2 m\ /dt\ 2 
&~r) \fo) 


We can immediately integrate (69.2) and (69.3) and the 
results are 



where h and A: are constants. On eliminating i and a 
from (69,4) and (69,5) we obtain 

1 /dr\ 3 1 / 2 m\ c 2 k 2 1 / 2m\ 

— ^ \ty) — r* \ l ~d*r) + T = F*\ 1_ 

Now substitute r = Iju and differentiate the equation so 
obtained with respect to y> and the result is 


(69.6) 


d 2 u m 

__ + „ + 


3mu 2 


For the planets of our solar system, the term mjc^h 2 is 
much larger than 3 mu 2 fc 2 , But when we neglect this latter 
term, we obtain Newton’s equation for the motion of a 
planet. Thus the first approximation to the solution of 

(69,6) is the Newtonian solution u — { 14 -e?cos^)}, 

where e is the eccentricity of the elliptic orbit and £ is 
the longitude of perihelion, A second approximation 
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to the solution can then be obtained in the form 
m 

u ~ i 1 + ecm (V~S — A£)} 7 where Zlf — 

This means that the major axis of the elliptic orbit is 
slowly rotating about its focus (the sun). The increase 
in corresponding to a complete revolution ip = 2 jz is 
thus Qnibcjc^h 2 , For the planet Mercury the advance of 
perihelion is calculated from this to be 42.9 seconds of 
arc per century. This agrees well with the observational 
figures of 43.5 seconds of arc per century. 


§ 70, Einstein's universe 

Einstein was led by cosmological considerations to con¬ 
sider the universe with the metric 


(70.1) dcr 3 =-—(I— r 2 j3l 2 )^dr 2 — r 2 dd 2 — r 2 sin* Odip 2 +c*di* t 


where ^ is a constant. This metric is spherically symme¬ 
trical with e ~ A “ (1 — f 2 /^ 2 ) and v = 0, The Christoffcl 
symbols of the second kind are readily obtained from 
(07,2), 

Let us investigate the path of a ray of light in Ein¬ 
stein’s universe. The path must be a null-geodesic and 
so its equations are given by three of the four equations 

(26,4) taken together with g u — -— = 0. That is. we 

have du du 




(7o - 2) S+ 7 ££“ sinecose 

d 2 w ( 2 dw dr dip dO 

d^ + ~rd^dH + COt0 dHd^ = 0, 

dH 

-(1 -r*/« , )-(*)*.- r .(®)‘-, !si „.s(|)‘ +c .(|) , _o, 
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where u is some parameter. Again, following the argument 
of section 69, the equation (70.2) tells us that we can 
take 0 to have the permanent value rc/2. With this choice, 
the remaining equations reduce to 


(70.3) 


d*y) 2 dip dr & % t 

du* r du du * du 3 


(70,4) — (I— 


■2\-l 


I^Y—rd^Y 

\du} \du) 




We integrate (70,3) and obtain 

dt 


(70.5) 


r 2 ?! = h — = k 

du n> du *' 



where k and k are constants, and then eliminate i and u 
from these equations and (70,4)* The result is 



The solution of this equation is 


II c* &* 

(™.6) fa - ^2 cos2 (V “ f) + ~tfT sina iV ~ * > 

where £ is a constant. We immediately sec that r regains 
its initial value when ip is increased by n and that r is 
never infinite for any value of ip , Thus all the nulb 
geodesics of Einstein’s universe, that is the light rays, 

dt k 

arc closed curves, W 7 e see from (70.5) that = “ r 2 , 

v ' dip h 

Hence the time taken for a light ray to make a complete 
circuit is given by 

h ( 1 r 2 k 2 I” 1 

(if.” 8 ’(<*-{> +-51-■>» , (v-f)) <¥ 
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Owing to the periodicity of ip we have 

k C 2it f 1 c 2 k 2 l -i 

T = irJo eos2 ^ + sh H dy> 

4* M C S fc* . 1-1 

= tJ 0 [w* cos v + ~w 8111 v j * 

and on carrying out the integration we find thatT=2 tt^/c. 

Ex* l. Show that the Einstein universe is neither an Einstein 
space nor a space of constant curvature* 

Ex* 2* Prove that the curvature invariant of Einstein^ universe 
is It = 6 /^s. 


§ 71. De Sitter’s universe 
Other cosmological considerations suggested to De 
Sitter that the universe could be described by the metric 

{71.1) do 2 = — (1 — r 2 /^ a )-i dr 2 — rW—r* sin* 6 dip* 

+ c*( I— r 2 J0l 2 )dt % . 

This metric has also spherical symmetry but the con¬ 
stant 0fc has not the same value as the corresponding 
constant of the Einstein universe. 

The paths of light rays are the null-geodesics given by 
the equations 


(71.2) 

(71.3) 


d?Q 2 dr dB . A 

xt H-r - 3 - sm A cos 0 

du* r du au 


d 2 tp 2 dr dip d6 dip 

77 +- r ~ + 2 cot 0 — ~ = 0, 

r du du du du 




du* 5 ?* 1 >' 




<71 - 4) (l)'- r '(S)'-’ Jsi " , 9 (S) ! 

+ «*(!-,■/«■) (I)’ - 0 . 
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Again we can choose 0 = ji/2 permanently, and on in¬ 
tegration the equations (71.2) and (71.3) become 

S-*(> 

We now eliminate 2 and w from these equations and 

(71.4) . The result is 

where a ! = c l k 2 jh 2 + 1 l&l 2 . This equation can be im¬ 
mediately integrated and yields 

(71.5) 1 =acos(y> — #), 

where f is a constant. These trajectories correspond to 
straight lines and are not closed, since r becomes infinite 
when ip — f = jr/2* 

Ex, Show tfiat the De Sitter universe is an Einstein space with 
constant curvature 12/#*, 

Solution 

p. 114* Ex, e -A — 1 + ar 4 -f b/r; v *= k — h where a, b and k 
are constants. 
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